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ELECTRON  DISTRIBUTION  FUNCTIONS  AND  MOBILITY  IN  III-V  SEMICONDUCTOR 
COMPOUNDS  FOR  HIGH  AND  LOV  ELECTRIC  FIELDS 

1.  INTRODUCTION 

Since  the  1950s  when  the  first  paper  was  published  about  the  semi¬ 
conducting  properties  of  III-V  compounds,  this  class  of  compounds  has 
been  used  In  various  electronic  devices  such  as  infra-red  detectors, 
semiconductor  lasers.  Gunn  diodes,  and  light  emitting  diodes.  The 
construction  of  these  devices  is  based  on  various  transport  properties. 
One  of  the  most  important  of  these  is  the  drift  mobility,  which  will  be 
discussed  in  this  report. 

Determination  of  the  drift  mobility  of  III-V  semiconductors  is 
useful  in  the  development  of  high  speed  devices.  In  this  report  the 
electron  drift  mobility  is  found  here  as  a  moment  integral  over  the 
momentum  spaoe  of  the  eleotron  distribution  functions.  The  distribution 
function  la  calculated  from  the  steady  state  Boltzmann  equation  whioh 
balances  the  frictional  forces  due  to  scattering  of  the  oarriers  and  the 
foroe  on  the  oarriers  due  to  the  applied  electrio  field.  The  properties 
of  the  eleotron  distribution  functions  are  of  interest  itself  as  well  as 
a  means  for  oaloulatlng  the  drift  mobility. 


In  the  III-V  semiconductor  compounds  polar-optical  scattering  is 
the  dominant  form  of  scattering  mechanism  over  a  considerable  range  of 
temperatureU-41 .  Since  polar  scattering  is  not  elastic,  a  relaxation 
time  cannot  be  easily  written.  A  relaxation  time  exists  only  when  the 
carriers  have  energy  which  is  much  larger  or  much  smaller  than  the  polar 
optical  phonon  energy[5,61.  Therefore  the  relaxation  time  approach  is 
not  suitable  for  solving  the  Boltzmann  equation  in  polar  compounds.  A 
few  methods  have  been  suggested.  The  earliest  of  these  is  the 
Variational  Method  [1,7,81.  Kohler[8l  introduced  this  method  in  1948; 
it  was  applied  to  polar  mode  scattering  by  Howarth  and  Sondheimer  [7] 
for  a  spherical,  parabolic  energy  band.  Ehrenreich[3]  later  used  the 
method  for  a  non-parabolic  Kane  [9]  band  where  he  combined  the  polar 
scattering  with  acoustic  and  electron-hole  scattering.  In  many  cases 
the  method  was  and  still  is  mathematically  complicated  when  results  of 
sufficient  accuracy  are  required.  This  is  especially  true  when  non- 
parabolicity  of  the  bands  has  to  be  included.  A  more  attractive  method, 
the  Drifted  Maxwell-Boltzmann  Distribution  Method  [10,111,  was  widely 
used.  To  some  extend  it  is  still  used.  This  method  assumes  that  the 
carrier  distribution  function  can  be  approximated  by  a  drifted 
Maxwellian  function.  However  this,  method  is  not  applicable  to  some 
III-V  semiconductor  compounds  at  high  electric  fields  [121. 

At  present  two  of  the  more  popular  methods  are  the  Iteration  Method 
[13-151  and  the  Monte  Carlo  Method  [17-211.  In  the  Iteration  Method  the 
Boltzmann  equation  is  reduced  to  an  infinite  set  of  algebraic  equations 
which  are  truncated  at  some  point  and  solved  iteratively.  This  method 
has  the  advantage  of  self-correction  of  the  error  through  iteration. 


Very  often  the  number  of  calculations  required  is  smaller  than  in  the 
Variational  Method.  Rode! IS. Id]  made  a  series  of  calculations  using 
this  method  and  was  able  to  show  that  its  results  agree  very  well  with 
experimental  data  for  a  number  of  III-V  compounds  such  as  GaAs.  GaSb. 
InAs  and  InSb  [15], 

Vlth  the  availability  of  better  and  faster  computers,  the  Monte 
Carlo  Method  [21]  has  become  more  popular.  In  this  method  the  Boltzmann 
equation  is  not  direotly  solved.  The  transport  coefficients  are 
evaluated  by  the  simulation  of  eleotron  trajectories  in  a  computer.  The 
solution  to  the  transport  equation  can  be  approximated  by  estimating  the 
characteristics  of  the  oarrler  population  as  statistical  averages  over  a 
large  number  of  collisions.  This  method  is  especially  useful  in  those 
oases  where  it  is  difficult  to  solve  the  transport  equation  by  the  other 
methods  described.  The  main  drawbacks  of  this  method  are  the  large 
amount  of  computer  time  required  for  the  statistical  averaging  and  the 
ad  hoo  assumption  about  the  collision  prooess. 

In  this  report  the  oaloulations  of  the  electron  distribution 
function  and  drift  mobility  are  made  with  some  assumptions.  These  are 
that  there  is  an  applied  DC  electric  field  and  that  the  carriers 
(electrons)  are  in  a  steady  state.  Under  these  conditions  seml-olassl- 
oal  statistics  are  often  used.  The  Boltzmann  equation  oan  be  used  if 
the  device  is  not  too  small  and  if  the  temperature  is  not  too  low.  For 
very  low  temperature  or  for  very  small  devioes,  quantum  effeots  have  to 
be  considered.  Calculations  will  be  made  for  low  electric  field  and 
high  electric  field  cases.  For  the  low  elec trio  field  case  two 
additional  methods  will  be  introduced.  These  two  methods  are  extensions 


-4- 


-•  _  * 


of  the  Iteration  Method  used  by  Fortini,  Dlquet  and  Lugand  [14].  The 
Iteration  Method  will  be  presented  briefly  to  verify  their  results  and 
to  use  them  as  a  basis  for  comparision.  It  Is  part  of  the  aim  of  this 
report  to  arrive  at  a  method  vhloh  is  less  tedious  than  the  Iteration 
Method  and  the  Monte  Carlo  Method.  The  Analytical  Method  is  especially 
of  interest  here  because  it  allows  the  expression  of  the  asymmetric  part 
of  the  electron  distribution  funotion  in  a  simple  and  compact  form.  The 
analytic  form  gives  a  better  picture  of  the  different  parameter 
dependences,  such  as  temperature  dependence,  of  the  distribution 
function. 

For  high  electric  fields,  a  simple  analytical  expression  for  the 
electron  distribution  function  has  been  derived  by  Stratton  [5].  and 
Conwell  and  Vassell  [6].  However,  one  of  the  drawbacks  of  their  analy¬ 
tical  result  is  that  it  fails  at  low  energy.  In  this  report  a  Pertur¬ 
bation  Method  is  suggested,  which  extends  their  result  to  include  the 
low  energy  region.  The  effect  of  non-parabolic ity  of  the  central  valley 
of  the  conduction  band  and  p-wave  admixture  on  the  carrier  distribution 
function  at  high  electric  fields  is  also  examined.  This  is  of  interest 
because  it  has  been  claimed  that  non-parabolicity  alone  can  cause 
negative  differential  mobility  within  the  central  valley  for  high 
electric  fields  [22]. 

1.1  BOLTZMANN  EQUATION  FOR  ELECTRONS  IN  DC  FIELDS. 

In  this  section  the  steady-state  Boltzmann  equation  for  electrons 
in  III-V  semloonduotor  compounds  with  a  non-parabolio  energy  band  and  p- 
wave  scattering  included  is  presented.  For  a  steady  state  system  the 
Boltzmann  equation  consists  of  two  terms,  a  collision  term  and  an 
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electric  field  term.  The  collision  term  represents  the  rate  of  change 
of  the  carrier  distribution  funotion  due  to  the  collisions  and  the 
electric  field  term  is  the  rate  of  change  of  the  carrier  distribution 
function  oaused  by  the  force  on  the  carriers  due  to  the  applied  electrio 
field  2. 

(k)e*(k)c-»- 

The  electric  field  term  and  the  collision  term  ar»  indicated  by  the 
subscripts  C  and  C  respectively.  In  this  report  carriers  are  assumed  to 
be  electrons. 

The  collision  Integral  consists  of  four  terms  that  represent  the 
scattering  in  momentum  space  of  electrons  in  and  out  of  the  state  with 
wave  vector  £.  The  terms  are  the  scattering  in  and  out  of  state  £  with 
emits  ion  or  absorption  of  a  phonon. 

[dt]  |N> l2LClc,ic (tc'> 

+1  (K-l \ H»  I H)  1 2L flc.lc ' ) S ,-Dw) f  (Ic ' ) 
-l(N+llH'lM)l2L(Ic.lc’)5(|K,-CR+H«)f(£) 

-I  (N-l  I H  •  1 N)  1 2L  (£ M ' )  6  (Eg  ,-gg-Du)  f  (K)J 

(1.2) 

D  is  the  Planok's  constant  divided  by  2n  and  u  is  the  frequency  of  the 
longitudinal  optloal  phonon.  E  ^  is  the  energy  of  the  carrier  asso> 
dated  with  the  state  ic  and  the  difference  in  energy  between  ic  and  k’  is 
Dm.  In  (1.2)  the  terms  represent  scattering  of  carriers  into  k  from  k’ 
by  entoion  of  Dm,  into  k  from  k*  by  absorption  of  Du,  out  of  k  by 
emission  of  Dm  and  out  of  k  by  absorption  of  Du,  respectively.  L(k,k*) 


The  term  (N+llH'lN)  in  (1.2)  Is  a  matrix  element  describing  the 
scattering  events  for  the  transition  between  the  states  k  and  ic’  brought 
about  by  a  perturbing  potential  H*.  For  high  mobility  materials  such  as 
those  considered  here,  the  coupling  between  the  electron  and  phonon  is 
weak.  Perturbation  theory  may  be  used  to  desoribe  the  interaction.  The 
matrix  element  for  this  oase  has  been  derived  by  Frolich  [24],  Callen 
[23]  and  Ehrenrelch  [1]  using  perturbation  theory.  The  matrix  element 
for  a  wave  function  with  s-symoetry  can  be  written  as  [26] 
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I  /N+^  +  AN\ 

m*Vlic-?c  •  I  »\  2  2/ 


(1.7) 


where  H  Is  the  number  of  phonons  with  wave  vector  k'-k  and  AN  is  the 
change  In  N  in  the  process  described  by  the  matrix  element.  V  is  the 
volume  of  the  crystal  and  eEQ  is  given  by  (26] 


eE.  -  a  e^hw 
V 


(Vi) 


(1.8) 


with  e  being  the  charge  of  an  electron  and  tm  and  \  are  the  dielectric 
constants  in  cgs  units  for  infinite  and  zero  frequencies  respectively. 


To  reduce  equation  (1.2)  to  a  more  useful  fora,  the  summation  over 
ic'  Is  transformed  to  an  integration  in  spherical  coordinates  with  k 
along  the  z-axis.  For  a  weak  electric  field  the  electron  distribution 
Amotion  is  nearly  spherically  symmetric ( 27 ] .  This  electron  distribu¬ 
tion  may  be  expanded  into  a  series  of  Legendre  polynomials  that  is  made 
up  of  a  dominant  term  and  a  series  of  small  perturbed  terms.  It  has 
been  shown  for  the  case  considered  here  that  retaining  the  dominant  term 
and  the  first  asymmetric  term  gives  a  good  approximation  to  the  electron 
distribution  function[28] .  The  distribution  function  in  the  collision 
integral  can  therefore  be  approximated  by  the  first  two  Legendre 
polynomials 


->•:  v 


f(C)  -  fQ(£)  *  kggte) 


(1.9) 


where  f0(£)  and  kgg(£)  are  the  symmetric  and  perturbed  parts  of  the 
distribution  Amotion  respectively.  From  now  on  the  subscripts  ic  and  k' 
will  be  dropped.  The  notations  £  and  £'  will  be  used  instead.  k£  is 
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the  oomponent  of  k  in  the  direction  of  the  eleotrio  field  £ .  When 
equations  (1.3)-(1.9)  are  substituted  into  the  Integral  fora  of  (1.2), 
it  ean  be  integrated  as  shown  in  Appendix  1.  The  result  of  the 
Integration  can  be  separated  into  a  symmetrical  part  and  an  asymmetri¬ 
cal  part.  The  two  parts  of  the  rate  of  change  of  the  electron 
distribution  due  to  collisions  are  shown  below: 

|df ^  w  eEQM  j(f0(t')e°  -f0(t))|A'2ln|^7^|-A'co'-aa'ce'j  dQ 
0  sya  C’-t+hu 

+(f0(C')-e°  f0(t))^A*2ln|g~|-Abo'-aa’ce'jdn  J 
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(1.10) 
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where 


A*  «  aa*  +  (k*)2  +  k2cc* 
2kk» 


(1.11) 

(1.12) 


and 

a  -  WkgT  ;  (1.13) 

kB  is  the  Boltzmann  constant  and  T  is  the  electron  temperature.  In 
(1.10)  and  (1.11)  Q  is  the  energy  function  defined  in  terms  of  the 


effective  mass. 

0  -  h2k2/2m*  -  (  1  +  at  +  HI2  +  ...  )  (1.14) 


The  coefficients  a,  p . 


see 


are  constants  that  can  be  calculated  for  each 


compound .  The  details  are  disoussed  in  referenoe  [6]. 

The  Boltzmann  collision  terms  (1.10)  and  (1.11)  contain  both  a  non¬ 
parabolic  energy  band  for  the  central  valley  of  the  conduction  band  and 
a -wave  and  p-wave  scattering.  Many  earlier  treatments  included  non- 
parabolio  energy  bands  with  s-wave  scattering.  However,  this 
approximation  is  inconsistent  because  a  non-parabolio  band  exists  only 
If  p-wave  and  s-wave  matrix  elements  are  present;  a  parabolic  energy 
band  is  consistent  with  only  s-wave  matrix  elements. 

The  other  term  in  (1.1),  (df/dt)^,  which  is  the  rate  of  change  of 
the  electron  distribution  funotlon  due  to  the  electric  field,  has  been 
derived  by  Conwell  and  Vasselfdl  and  is  given  by 


In  section  2,  the  electron  distribution  function  is  found  from  the 


steady-state  Boltzmann  equation  for  low  electric  fields  by  the  Iteration 
Method . 


For  low  electric  fields,  there  is  very  little  mixing  of  the  s-wave 
with  the  p-wave.  Therefore  c  and  o'  in  (1.10)  and  (1.11)  can  be  approx* 
i mated  by  zeros,  k’  is  reduced  to  unity  in  the  same  equations.  It  is 
assumed  in  this  section  that  the  energy  band  is  spherical.  Therefore  0 
is  equal  to  £  and  dQ  /d£  is  equal  to  unity  in  (1.10)-(1.12) .  Under  this 
condition  the  symmetrical  part  of  the  Legendre  expansion  of  the  electron 
distribution  function,  fQ,  oan  be  approximated  by  a  thermal  equilibrium 
Max we 11-Bolt zraann  distribution. 
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(2.1) 


where  is  the  effective  density  of  state  and  nQ  is  the  density  of 
electrons  in  the  conduction  band  valid  for  low  electric  fields.  Clas¬ 
sical  statistics  for  the  electrons  are  used  here.  Under  certain  condi¬ 
tions  Ferml-Dlrao  statistics  must  be  used. 

With  (2.1),  (1.12),  which  is  the  asymmetrical  part  of  the  field 
equation,  can  be  equated  with  the  asymmetrical  part  of  the  collision 
term  according  to  the  steady-state  Boltzmann  equation  to  give 
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where  H($-e  )  is  the  unit-step  function,  $  is  the  normalised  form  of  the 
energy,  E/lc^T.  This  equation  is  the  same  as  equation  (S)  of  the  paper 
of  Fortini  et  al  (141  except  there  is  a  missprint  in  their  paper. 

The  square-root  sign  in  the  denominator  is  over  (§  -e  )  instead  of  the 
product  of  ?(?-«).  If  a  new  function,  y(5,e),  is  defined  as  shown  in 
equation  (2.3),  equation  (2.2)  can  be  rewritten  in  the  form  of.  equation 
(6)  of  the  paper  by  Fortini  et  al  Cl4}.  This  is  shown  in  equation  (2.4) 
below. 
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where  the  ooefiolents  A($),  BO|)»  C({),  and  D($)  are  defined  as 
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2.1  ITERATION  EQUATIONS  FOR  THE  ASYMMETRICAL  DISTRIBUTION  FUNCTION 
One  of  the  methods  for  solving  the  difference  equation  (2.4) 
Involves  iteration.  This  has  been  done  by  Fortini,  Diquet  and  Lugand 
(141,  and  Grigorev,  Dykman  and  Tomohuk  (13].  In  this  section  the 
method  of  Fortini  et  al  (14]  will  be  described  briefly.  Although  the 
aim  of  this  report  is  to  find  approximate  analytical  expressions  for  the 
electron  distribution  funotion,  the  Iteration  Method  is  solved  in  this 
section  in  order  to  verify  the  results  of  Fortini  et  al.  The  results 
then  can  be  used  as  a  reference  by  which  the  resuits  from  the  other 
methods  may  be  compared. 

A  new  variable  {  is  defined  such  that  $  »  5  +n<r  where  the  range 
of  ?  is  from  zero  to  a  .  Therefore,  for  a  given  value  of  a  the  distri¬ 
bution  funotion  consists  of  a  set  of  functions  y (?,«  ),  y(£+<j  ,a  ),  .... 
y ( 5+nc  ,a  ),  Each  of  these  y(S+ne  ,o  )  is  defined  in  the  interval  no 
to  (n+l)o  .  For  convenience,  y($+ne  ,e  )  will  be  denoted  by  yn<  The 
funotion  y($,e  )  may  be  found  for  the  range  of  (  needed  by  combining  the 
set  of  yQ  for  n  equals  to  0,  1,  2,  ...etc.  In  terms  of  yR  equation 
(2.4)  can  be  written  as  an  infinite  set  of  coupled  difference  equations: 

Vo  *  Vl  ■  Do 
*1Y0  +  B1Y1  +  C1y2  "  D1 
Vl  +  B2y2  +  C2y3  “  °2 

(2.9) 

*n-lY n-2  +  Bn-lYn-l  +  Cn-lYn  “  Dn-1 
*nYn-l  +  ®nYn  +  CnYn+l  "  Dn 

e 

e 

where  the  notation  for  A,  B,  C  and  D  is  defined  by 


A  finite  set  of  equations  may  be  ahosen  from  (2.9)  to  calaulate 
y($»o  ).  Since  there  Is  always  one  more  unknown  than  algebriae  equa¬ 
tions,  an  approximation  is  needed.  One  of  the  most  obvious  approxima¬ 
tions  is  the  assumption  that  cn7n+1  is  small  in  the  last  equation  of  the 
chosen  set  and  oan  be  dropped.  This  assumption  allows  us  to  truncate 
the  Infinite  set  of  coupled  difference  equations..  The  result  is  a  set 
of  n+1  equations  with  n+1  unknowns: 


These  equations  can  be  solved  simultaneously  for  the  unknown  yQ,  yJf 
r2****rn*  For  a  given  a  ,  y ($,e  )  oan  be  calculated  once  the  set  of  yQ 

has  been  found.  This  procedure  is  repeated  with  more  equations  (larger 
n)  until  the  change  in  the  yR  of  interest  is  small.  This  iteration 
procedure  is  expected  to  produoe  yn  to  any  desired  degree  of  aoouraoy. 
However,  the  oaloulation  becomes  very  tedious  as  n  increases. 


2.2  RESULTS  POR  y($.<r)  FROM  THE  ITERATION  METHOD 

With  the  term  CnTn+1  in  the  equation  (2.8)  dropped,  the  resulting 


set  of  n+1  equations  with  n+1  unknowns  may  be  solved  on  a  computer. 
Portinl  et  al  [14]  performed  the  iteration  on  an  IBM  computer  up  to  an 
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order  n  where  rn+1  differs  from  Yn  by  less  than  5x10  .  y(?.o  )  was 

obtained  in  the  range  0<£<(n+l)®  from  the  set  (Yq.Yj.Y2* • • «Yn)  • 

However  no  values  for  n  were  given  to  indicate  the  rate  of  convergence 
for  the  Iteration  process. 

In  this  report  the  iteration  is  done  on  an  Apple  II  computer  and 
the  results  of  y(Z.a)  are  plotted  against  {  in  figure  1.  The  plots  are 
for  values  of  «r  equal  1.  2,  3,  4.  and  5.45.  These  results  are  almost 
identical  to  those  given  by  the  Fortini  et  al  [141  paper.  In  that  paper 
only  plots  for  a  equal  1,  2,  3  and  4  are  given.  The  additional  plot  of 
o»5.45  in  this  report  corresponds  to  the  temperature  of  liquid  Nitrogen, 
77®K.  This  value  is  chosen  because  of  the  interest  in  devices  working  at 
that  temperature.  The  temperature  equivalent  to  the  optical  phonon 
energy,  Hu.  is  assumed  to  be  418°*  as  compared  to  416°K  used  in  the 
reference  [141.  Room  temperature  will  correspond  to  a  a  of  1.393  and 
is  expected  to  have  a  plot  of  y($.®  )  lying  between  a  -1  and  a  -2.  At 
values  of  {  where  a  is  an  integer,  there  are  sharp  breaks  because  of 
the  emission  of  a  longitudinal  optical  phonon.  These  sharp  breaks  are 
not  discontinuities.  They  appear  discontinuous  because  of  the  scale 
used  in  plotting  y(£»®  ).  An  expansion  of  the  scale  around  these  'break 
points'  will  show  that  they  are  continuous. 
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dotted  line  is  for  e- 2.  the  dotted  line  is  for  <r»3,  the  large  dashed 
line  Is  for  o-4  and  the  small  dashed  line  is  for  o«5.45. 


To  get  an  idea  of  the  rate  of  convergence  of  the  iteration 
process,  values  of  yQ  at  £*0.7  are  calculated  for  increasing  number  of 
equations,  n.  The  result  are  tabulated  in  table  1.  The  value  g*0.7  is 
ohosen  because  it  is  near  the  peak  of  the  distribution  funotion  and  the 
absolute  difference  is  expeoted  to  be  greatest  there.  Prom  the  table  it 
is  apparent  that  the  convergence  for  the  iteration  is  slow  for  cr  *1,  but 
as  a  increases  the  rate  of  convergence  increases  sharply.  For  a  *1, 
five  equations  (n*4)  are  needed  for  the  ohange  of  0.3  percent  in  y  ,  but 
for  or  *4  the  same  number  of  equation  will  give  a  difference  of  less  than 
5xlO-5  percent  in  y_. 

In  table  2,  values  of  yR  at  £*0.7  for  various  <i  are  tabulated. 

This  table  gives  a  comparison  of  the  magnitude  for  y  of  various  a  . 

n 

This  allows  the  estimation  of  the  errors  resulting  from  the  truncation 
of  the  infinite  set  of  equations  into  a  set  of  n  equations.  From  these 
two  tables,  it  is  apparent  that  for  large  values  of  a  ,  yQand  y.  alone 
will  provide  a  good  approximation  for  y(£,o  ). 


2 


3 


5.45 


DBS 

i 

- 

0.39263 

0.43313 

0.47068 

0.52592 

2 

- 

0.40291 

0.43603 

0.47142 

0.52603 

3 

0.36417 

0.40403 

0.43613 

0.47143 

0.52603 

4 

0.36593 

0.40415 

0.43614 

0.47143 

0.52603 

5 

0.36638 

0.40416 

0.43614 

6 

0.36649 

7 

0.36652 

3 

0.36653 

Table  1.  Values  of  different  at  $-0.7  when  different 
number  of  equations  are  used.  The  number  of  equations  used 
in  each  row  is  n+1. 
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The  Boltzmann  equation  for  low  eleotrlo  fields  la  solved  by  the 
Perturbation  Method  In  section  3.  The  Perturbation  Method  is  introduced 
as  a  first  step  in  finding  an  approximate  analytical  form  for  the  elec¬ 
tron  distribution  functions. 


In  order  that  a  perturbation  expansion  can  be  done,  a  small  expan¬ 
sion  parameter  is  needed .  e~°  can  be  used  since  it  is  smaller  than 

unity  for  all  positive  a  .  For  most  of  the  temperature  range  of 
Interest,  e~9  is  small  enough  that  rapid  convergence  can  be  expected. 

In  this  report  e“®  will  be  represented  by  5.  It  is  assumed  that  each  y. 


can  be  expanded  into  a  power  series  in 
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(3.1) 


The  superscript  of  y~  indicates  the  order  of  the  term  in  each  expansion 

n 

and  not  a  power.  The  coefficients  A  ,  B  .  C  and  D„  are  replaced  by  A®  , 

nan  n  n 

B°+B-1/&,  C"1/*  and  Dn5n  respectively.  In  this  form  the  coefficients  are 
n  n  n  n 

shown  below: 
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When  (3.1) -(3.6)  are  substituted  Into  (2.11),  the  result  is  the  set  of 
equations  shown  In  (3.7). 
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Relations  between  the  terns  In  the  expansions  may  be  obtained  by 

n 

equating  the  coefficients  for  various  powers  of  &.  By  equating  the 
coeflclents  of  ft”1,  a  set  of  equations  relating  y®  for  different  n  is 
obtained: 
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Prom  the  first  equation  in  (3.8)  ,  y®  is  zero,  y®  «0  implies  all  y® 
with  n  greater  than  zero  are  zeros.  This  oan  be  seen  by  substituting  y^ 

into  the  seoond  equation,  finding  y®-0  and  then  substituting  y®  into  the 
next  equation,  etc.... 
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for  n>0 


(3.9) 
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A  set  of  equations  relating  and  y*  can  be  obtained  by  equating  the 
coefficients  of  81.  This  set  of  equation,  which  is  shown  in  (3,10),  can 
be  reduced  to  the  two  equations  in  (3.11)  by  using  the  results  for  y® 
obtained.  The  two  equations  in  (3.11)  oan  be  solved  simultaneously  for 
rj}  and  yj: 

*0  0  +  C”1*1  -D° 

Vo  0  Y1  D0 

*0  0  +  a  a  n-K1  m  ft 

Vo  +  iYl  +  B1  Y1  +  C1  y2  0 


(3.10) 
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The  procedure  for  obtaining  expressions  for  other  y"  is  very  simi¬ 
lar  to  that  use  for  y^  and  y*  .  In  general,  by  equating  the  coeffi¬ 
cients  of  8°,  a  set  of  equations  relating  y®-1  and  y®  is  obtained.  This 
set  of  equations  can  be  reduced  to  m+2  equations  by  using  the  results 
that  have  been  found  from  solving  the  equations  obtained  by  equating  the 
coefficients  of  5n_1  .  The  set  of  m+2  equation  is  shown  in  (3,12). 

❖Z  *  W  -  ° 

‘ft  *  *ft  <rr‘  *  ■  • 


,0  m-1  .  o0  m-1  .  _-l  m  ,  _-l  m  A 
A  ,y  a+B  ,y  ,  +  B  ,y  «+C  ,y  *0 
m-lrm-2  m-1 1 m-1  m-lfra-l  o-lrm 


.0  m  ^  _0  m  ,n-l  m+1  .  --1  m+1  -m 

Ay  «  +  B  y  +By  +C  y  ,  =  D 

m  m-1  m  m  m  'ra  m  m+1  m 


(3.12) 


■  -(v£K’rl<Ci)-  »«•  <3-io> 


Unlike  the  Iteration  Method,  the  aocuraoy  of  any  rn  calculated  by 
the  Perturbation  Method  does  not  depend  on  the  truncation  process.  The 
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accuracy  of  the  Perturbation  Method  depends  on  the  number  of  terms 

summed  for  each  yq.  When  equating  the  coefficients  for  a  given  Rm  in 

(3.7),  so  long  as  the  number  of  equations  n  is  picked  such  that  n>m+2, 

the  set  of  equations  will  reduce  to  just  m+2  equations.  With  these  m+2 

equations,  the  m+2  unknowns  in  the  expansions  can  be  calculated 

n 

exactly. 


Theoretically  the  Perturbation  Method  allows  us  to  obtain  a  given 
yr  to  any  degree  of  accuracy  by  using  (3.13)-(3.20) .  It  is  apparent 
that  to  calculate  a  oertain  term,  a  few  terms  of  the  preceding  order 
have  to  be  used.  However  each  term  needs  at  most  three  substitutions. 
Table  3  shows  the  terms  collected  after  solving  the  set  of  equations 
obtained  from  equating  each  power  of  5  in  (3.7).  The  terms  in  each 
column  are  solved  simultaneously  by  using  the  terms  in  the  column  to  the 
left. 


isVW, 


•  .*■ 


Table  3.  The  different  rn  in  each  oolumn  are  the  terms  that 
oan  be  calculated  from  the  set  of  equations 
obtained  by  equating  the  coefficients  for  the  S 
shown  above  each  column. 


3.1  y($,<i  )  FROM  THE  PERTURBATION  METHOD. 

Plots  of  y(?.«  )  versus  $  obtained  by  the  Perturbation  Method  are 
shown  in  figure  2.  The  plots  are  for  a  -1.  2,  3,  4  and  5.4S.  These 
plots  are  almost  identical  to  those  obtained  by  using  the  Iteration 
Method  except  for  (  near  zero.  Different  number  of  terms  are  summed  for 
eaoh  rQ.  for  a  -1,  the  number  of  terms  summed  for  y0  is  nine,  and  for 
other  a  ,  the  number  of  terms  summed  is  10-n.  For  9  *4  .  three  and  two 
terms  are  summed  for  yq  and  Yj  respectively.  When  9  is  large,  the 
convergence  is  very  rapid  and  only  a  few  terms  in  each  Yn  expansion  are 
required  to  provide  a  good  approximation. 

Comparisons  between  the  sizes  of  some  of  the  terms  are  shown  in 

table  4  and  3.  The  comparisons  are  done  at  1*  >0.7.  In  table  4  one  of 

the  observable  trends  is  that  for  n>0,  the  dominant  term  is  the  second 

term  in  each  expansion.  Even  though  this  is  a  little  unusual  in  a 

convergent  series,  it  is  not  unheard  of.  This  type  of  behaviour  can  be 

found  in  some  expansions  in  Bessel  functions  and  some  solutions  to  the 

Schrodlnger  equation.  It  is  apparent  that  the  first  three  terms  in  a 

given  y_  expansion  will  give  a  rather  good  approximation  for  y„. 
n  n 
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Table  4.  Values  of  the  different  terns  in  the  expansions.  The 
values  of  r  from  the  Iteration  Method  have  been  inoluded  in  the 
last  row  for  comparison. 
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Table  6  shows  the  size  of  the  errors  when  the  first  two  terms  and 

the  first  three  terms  are  used  to  approximate  v  .  The  results  from  the 

n 

Iteration  Method  of  section  2  have  been  used  as  standards.  The  compari¬ 
sons  are  made  for  a  -1  and  er  *4.  For  each  rn»  the  results  from  summing 
the  first  two  and  three  terms  are  divided  by  the  results  of  the  Iteration 
Method.  For  a  -1  the  error  from  using  three  terms  is  always  less  than 
one  percent.  The  error  from  using  two  terms  is  greater.  However  in 

most  cases  the  first  two  terms  alone  are  sufficient  to  represent  y  . 

n 

Sven  though  the  error  is  about  ten  percent  for  y  ^  for  the  case  when  a 

>1.  this  will  not  really  affect  the  mobility  calculations  very  much. 

This  is  because  the  contribution  to  the  mobility  by  each  y  becomes 

n 

smaller  as  a  gets  larger.  The  main  contributions  to  the  drift  mobility 
come  from  yQ  and  Tj. 

Since  the  leading  two  terms  are  enough  to  give  a  good  approximation 
for  a  given  y  ,  the  Perturbation  Method  is  simpler  than  the  Iteration 
Method . 


Ratio 


Table  6.  Ratio  of  the  y  from  the  Perturbation  Method  to  the  yn 
calculated  by  IterationVlethod.  The  ratios  tabulated  are  for 
a  *1  and  4,  at  5*0.7.  For  each  a  the  results  from  using  two  and 
three  terras  in  the  expansions  are  compared. 


3.2  FAILURE  OF  THE  PERTURBATION  METHOD 


For  g  near  zero,  the  Perturbation  Method  fails  to  approximate  y  . 

n 

The  magnitude*  of  the  yR  are  very  large  when  £  approaches  zero.  The 
failure  is  more  apparent  for  a  -1  but  is  also  true  for  the  other  a  . 

The  sign  of  the  yn  can  either  be  positive  or  negative,  and  can  be 
different  when  a  different  number  of  terms  Is  summed.  This  suggests 
that  the  expansions  for  yQ  are  divergent  at  ?  near  zero. 

For  £  small,  6  is  no  longer  small  compared  to  £.  The  Perturbation 
Method  fails  when  this  happens  because  terms  kept  for  the  lowest  order 
in  yR  are  smaller  than  terms  dropped.  This  result  can  be  most  easily 
seen  by  truncating  equation  (2.10)  to  two  equations.  The  results  for  yQ 

and  y4  from  solving  the  two  equations  are 


r0  “ 


B1D0  ‘  C0D1 

Vi  -  C0A1 


VrVo 

BqBj  -  kxC0 


(3.21) 


(3.22) 


When  both  6  and  £  are  small,  two  expansions  can  be  made  of  yQ  and  y^:  an 
expansion  in  6  and  an  expansion  in  ?.  If  5  and  (  are  very  small,  only 
the  leading  term  of  these  quantities  in  eaoh  expansion  has  to  be  kept. 
Since  two  expansions  in  a  small  parameter  are  made,  the  order  in  which 
the  expansions  are  performed  can  make  a  difference.  To  check  the  diffe¬ 
rence  an  expansion  first  in  5  and  then  one  in  £;  is  performed,  and  next 
the  order  in  whioh  the  expansions  are  done  is  reversed.  The  oheok  is 
done  on  yQ  and  y^.  To  obtain  the  leading  terms  in  5,  the  limit  as  6 
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approaches  zero  is  taken  for  y^  and  y^/R.  Por  the  leading  terns  in 
the  functions  that  approximate  yQ  and  yj  are  expanded  in  and  the 
leading  terns  are  extracted.  The  results  for  the  two  different  orders 
in  which  the  leading  terms  are  obtained  are  given  in  table  7. 


T0 

*1 

4- 

lioit5^rilmit^0(yn)] 

-/ST 

2 

gg 

Table  7.  The  limiting  forms  of  y-  and  y^  when  different 
order  in  limiting  procedures  are  used. 


Both  yQ  and  y1  approach  different  limiting  forms  for  the  two  limiting 
procedures.  However  for  yQ  the  error  is  not  so  significant  since  in 
both  limiting  procedures  yQ  approaches  zero.  The  difference  in  y^  is 
more  significant  since  the  first  limiting  procedure  leads  to  zero 
Instead  of  a  constant.  The  Perturbation  Method  gives  a  result  which 
corresponds  to  the  limiting  procedure  in  row  one  of  table  7.  For  ?  near 
zero  the  limiting  procedure  in  the  second  row  is  correct.  Therefore  the 
Perturbation  Method  fails  near  S?  equal  zero.  The  behaviour  of  other  y 

n 

can  be  examined  in  the  same  way.  One  would  find  that  the  Perturbation 
Method  fails  near  ?  equal  zero  as  for  yQ  and  y ^ . 

The  Perturbation  Method  fails  when  the  leading  terms  of  the  two 
limiting  procedures  are  of  the  same  order.  In  the  case  of  y^,  this 


occurs  when 


Since  8  is  a  decreasing  exponential  of  a  ,  as  a  increases,  the 
right  hand  side  of  (3.23)  decreases  rapidly.  Therefore  as  a  increases, 
the  point  where  the  perturbation  fails  moves  closer  to  g*0.  This 
explains  whv  the  failure  does  not  show  up  in  the  yR  plots  for  «  greater 
than  one.  If  the  horizontal  scale,  g,  is  expanded  near  g  «0.  the 
failure  is  expected  to  appear. 

Since  the  value  of  the  mobility  for  the  carriers  is  proportional  to 
the  area  under  the  y-curve,  the  Perturbation  Method  still  oan  give  a 
good  approximation.  This  result  follows  because  the  region  where  the 
perturbation  falls  badly  represents  only  a  small  fraction  of  g  between 
zero  and  a  .  The  small  section  oan  be  omitted  in  the  integration  of  y 

n 

without  seriously  affecting  the  values  of  the  mobility  calculated.  As 
long  as  the  integration  does  not  include  the  point  g  -0,  the  results  of 
the  mobility  should  be  good,  especially  when  a  is  large. 
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4  INTRODUCTION  TO  ANALYTICAL  METHOD 

Analytical  expressions  for  y($.<r  )  oan  be  derived  via  two  different 
approaches.  These  two  approaches  will  be  discussed  in  this  section. 

The  first  approach  is  to  derive  the  different  directly  from  the 
truncated  set  of  equations  in  (2.10).  Por  a  given  n,  this  set  of  equa¬ 
tions  oan  be  solved  sioultaneously  for  eaoh  v  in  terms  of  A  ,  B  ,  C  , 

r  nan 

and  Dn.  Efforts  oan  be  made  to  slmplifiy  the  result  of  each  yn  to  a 
simple  analytical  fora.  One  of  the  problems  of  this  approach  is  the 
complexity  involved  when  more  equations  in  (2.10)  are  used.  However, 
for  accurate  results  a  large  number  of  equations  have  to  be  used.  This 
will  lead  to  an  expression  for  ?n  whioh  is  very  complicated  and  cannot 
be  simplified  easily. 

The  second  approaoh  is  to  derive  the  analytical  formula  from  the 

result  of  the  perturbation  method  in  section  3.  In  the  perturbation 

method,  the  formula  for  any  y®  in  a  given  yR  expansion  has  been  derived 

and  is  given  by  one  of  the  formulas  in  (3.13)-(3.20) .  Each  of  these  y® 

n 

is  simplified  separately  and  they  are  combined  together  for  each  yn. 

This  method  is  expected  to  be  more  efficient  because  in  many  oases,  only 
the  first  few  y®  are  necessary  for  sufficient  accuracy.  Moreover,  it  is 
easier  to  simplify  eaoh  y®  than  the  yR  in  the  first  method. 


4.1  ANALYTICAL  FORMULA  FROM  THE  ITERATION  METHOD. 

In  this  section  the  method  of  deriving  analytical  formulas  for  y 

n 

by  solving  the  truncated  set  of  equations,  (2.10),  is  discussed.  The 
initial  step  is  to  decide  the  number  of  equations  to  be  used.  In  this 


method  the  aoouraey  of  the  result  for  the  eleotron  distribution  function 


depends  on  the  number  of  yR  used  and  the  acouraoy  of  eaoh  yr. 

Therefore,  as  many  equations  as  possible  should  be  included  in  solving 
for  the  different  y  ,  However,  the  number  of  equations  used  is  limited 
because  a  point  is  reached  when  the  results  for  these  yQ  found  by 
solving  the  set  of  equations  will  be  too  complicated  to  be  solved 
easily.  In  this  report  only  the  first  four  equations  will  be  used. 

If  (2.10)  is  truncated  at  n-3,  the  four  equations  are  shown  below: 

Vo  +  Vl  “  D0 
Vo  +  Vl  +  Clr2  *  D1 

Vl  +  B2Y2  +  C2Y3  *  °2  (4,1) 

A3y2  +  B3Y3  “  °3* 

These  equations  are  solved  simultaneously  and  the  expressions  for  y0# 
Yl’  y2*  and  are  shown  in  (4.2)-(4.S): 


can  be  grouped  In  terms  of  5.  Consider  the  expression  for  yQ  in  (4.2). 

The  A  B  .  C  and  D  are  replaced  by  A®  ,  B^+B^/S,  and  C^/S  and  o” 
n  n  n  n  n  n  n  n  n 

respectively.  These  factors  have  been  defined  in  (3.2)-(3.6) .  Since  6 

is  small  for  most  of  the  range  of  a  .  it  is  assumed  that  only  those 
—2  —3 

terms  with  S  and  6  dominate.  Other  terms  are  assumed  to  be 
relatively  small  and  dropped.  This  approximate  form  of  yQ  with  terms 
grouped  in  power  of  5  is  shown  in  (4,6). 

(4.6) 

Since  the  objective  here  is  to  arrive  at  a  simple  analytical  expression 

for  y0.  (4.6)  has  to  be  simplified  further.  A  new  variable  y ). 

is  introduced  for  the  convenience  of  the  simplification  process.  The 

Inverse  hyperbolic  sine  functions  in  A0,  B°,  B_1.  and  C-1  are  expanded 

nan  n 

into  Taylor  series  about  y-0.  Further  simplications  are  made  so  that 
they  can  be  written  in  polynomials  in  powers  of  y.  These  polynomials  are 
substituted  into  (4.6)  and  the  numerator  and  the  denominator  are 
simplified  separately.  A  relatively  simple  expression  for  yQ  can  be 
obtained  by  retaining  the  few  terms  with  low  order  in  y  and  5.  The 
final  form  of  yQ  is  shown  in  (4.7).  The  details  of  the  simplification 
process  are  shown  in  the  Appendix  3. 


-(5®  )1/2e"*  ry(1.01+0.4912y2)  +  1.56378' 


y(1.01-0,1260y  j  +  0.89048. 


(4.7) 


The  same  procedure  is  used  to  obtain  analytical  expressions  for  y^  and 
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2* 


-o  e 


0.67336  +  0.4488y' 
0.3366y  +  0.29686 


(4.8) 


y  m  __ 

I  ■)  mil  9  - 


3.37486  4-  4.0542y 


(4.9) 


Only  the  first  three  yQ  will  be  used  in  this  report. 

4.2  ANALYTICAL  PORMULA  FROM  THE  PERTURBATION  METHOD. 

In  this  section  the  analytical  foraula  for  rn  will  be  derived  from 
the  results  of  the  Perturbation  Method.  From  section  3,  it  is  apparent 
that  for  the  range  of  a  of  interest  here,  the  first  two  terms  or  three 
terms  in  each  rn  expansion  will  provide  a  very  good  approximation  for 
each  y  .  Therefore  only  the  analytical  expressions  for  the  first  few  rn 

need  to  be  found. 


The  first  analytical  expression  to  be  solved  is  the  leading  term  of 
the  yQ  expansion.  The  formula  for  y®  is  given  in  (3.13).  After  the 
substitution  of  A,  B,  C,  and  D  from  (2.5)-(2.8),  the  yjj  written  in  terms 
of  y  is  shown  below 


2sinh”1(y)5e”? 


[2sinh_1  (y)  ]  2-ri-(2y2+i)ainh‘'1  (y)  l2 


(4.10) 


/[y»(y»+l)J 


The  inverse  hyperbolio  sine  is  expanded  into  Taylor  series  about  y*>0, 
and  the  result  can  further  be  simplified  as  shown  in  the  Appendix  3. 

The  final  form  of  y^  written  as  a  polynomial  in  g/e  is  shown  in  (4.11) 


(4.11) 


In  this  oase  only  the  first  two  non-zero  tarns  have  been  shown.  With  y. 


found,  it  oan  be  used  to  find  yj  by  using  equation  (3.15).  The  A®  and  B”1 
in  that  equation  are  expanded  into  Taylor  series  about  x»0  as  before.  A 

polynomial  similar  to  (4.11)  is  obtained  for  y*  by  simplifying  the 

2  3  r 

resulting  expression.  This  prooess  oan  be  repeated  for  y*,  yjj,  . ..y 

12  3 

through  equations  (3.15W3.20) .  The  results  for  y*,  y*  and  r3  are 
shorn  below: 


fn 


rl 


-?e"^(l/3  -  0.05925?/®  +  ...  ) 
-?VC  (0.1544  -  0.01715?/®  +  ...  ) 

A* 

-?e"* (0.09505  -  0.00524?/®  +  ...), 


(4.12) 

(4.13) 

(4.14) 


In  order  to  obtain  y*  from  equation  (3.15).  the  results  for  y*  and 

Yj  nay  be  used.  As  in  the  previous  oases,  the  ooeffioients  a|,  C”1,  B® 

and  B"1  are  expanded  into  Taylor  series.  When  these  series  along  with 

(4.12)  and  (4.13)  are  substituted  into  (3.15),  the  resulting  expression 

oan  be  simplified  to  the  fora  shown  in  (4.15).  The  same  procedure  is 
2  3 

repeated  for  and  y2  by  using  equations  (3.18)  and  (3.19)  respectively 


1—2  M 

yj  -  -e  5(®  /3  +  0.1957?+  ...) 

y\  -  -®  e“*(1.345  +  0.2242?/®  +  ...) 

y\  -  -®  e“* (0.8724  +1.3052?/®  +  ...  ). 


(4.15) 

(4.15) 

(4.17) 


Before  higher  order  y*  are  oaloulated,  it  is  useful  to  see  how  eaoh  y* 

n  n 

will  oontribute  to  the  drift  mobility  calculation.  From  figure  1,  it  is 
apparent  that  yQ  and  y^  will  oontribute  most  signifloantly  to  the  drift 
mobility.  As  n  gets  larger,  the  contribution  from  a  given  yR  becomes 
smaller.  When  ®  is  greater  than  two,  the  contribution  is  mainly  from 
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yq  alone.  Therefore  It  is  not  necessary  to  oaloulate  rQ  to  more  than 
the  first  few  n.  For  eaoh  y  ,  it  can  be  seen  from  Table  6  that  if  only 
the  leading  two  terms  in  each  expansion  are  used,  the  approximations  are 
very  good.  In  order  to  get  simple  and  compact  analytical  formulas,  only 
the  two  leading  terms  are  used  for  eaoh  expansion  of  y^,  y2  and  y3. 
These  analytical  formula  are  shown  below: 

Yn  -  — (g/o  )1/2e~*/l  +  ll|\  -  e"?“<I(q/3  +  0.1967?)  (4.18) 

2  \  18«r/ 

Y.  -  -ge~*“a(. 3333-0. 0592g/q)-qe”5_2<,(l.  1346+0. 2242l/a) 

(4.19) 

y,  -  -^e”4“2o(0. 1644-0. 01716g/q)-q  e”*“3<T(0. 8724+1. 305g/o  ) 

(4.20) 

y3  -  -ge~*"3<T(0. 09605-0. 00S24g/q).  (4.21) 

4.3  y(g.q)  FROM  THE  ANALYTICAL  METHODS. 

The  resulting  yR  from  the  analytical  formulas  derived  from  the 
Iteration  Method  and  the  Perturbation  Method  are  plotted  in  figures  3 
and  4. 


In  figure  3,  the  yR  are  calculated  from  equation  (4.7)— (4.9) ,  and 

the  results  plotted  against  g  for  a  equal  to  1,  2,  3,  4  and  5.45.  There 

are  some  notioeble  differences  between  the  y_  in  figure  3  and  those  in 

n 

figures  1  and  2.  This  is  especially  true  for  smaller  a  .  The  analyti¬ 
cal  formulas  that  are  derived  from  (2.10)  oan  be  good  approximations 
when  a  is  large.  However  when  it  is  small,  the  approximation  is  not  as 
good. 


Methi 


The  problem  In  trying  to  derive  analytical  formulas  directly  from 
the  set  of  equations  in  (2.10)  is  the  complexity  involved  in  simplifying 
equations  (4.2)-(4.S).  As  dlsoussed  earlier  a  larger  number  of  equa¬ 
tions  in  (2.10)  are  needed  for  more  aocurate  y  ,  but  it  will  lead  to 
some  very  complicated  expressions  that  are  hard  to  simplify.  In  this 
ease  only  four  equations  have  been  used  and  yQ  is  very  complicated  as 
given  in  (4.2).  In  reducing  (4.2)  to  any  simple  analytical  form,  some 
of  the  terms  dropped  turn  out  to  be  significant.-  There  appears  to  be 
no  alternative  because  there  is  no  single  expansion  parameter  that  will 
enable  (4.2)  to  be  simplified  easily.  Therefore  simplicity  is  achieved 
at  the  expense  of  complexity. 

In  figure  4  the  result  of  the  yn  calculated  from  the  analytical 
formulas  derived  from  the  Perturbation  Method  are  plotted  against  5  for 
<r  equals  1.  2,  3.  4  and  5.45.  The  results  obtained  from  these  analyti¬ 
cal  formulas  are  very  close  to  those  obtained  in  figures  1  and  2.  The 
only  noticeble  difference  is  at  5-0  where  the  calculated  from  the 
analytical  formula  in  this  case  does  not  go  to  zero.  The  failure  of  the 
Perturbation  Method  as  discussed  in  section  3  is  the  root  of  this  prob¬ 
lem.  Since  the  drift  mobility  is  proportional  to  the  area  under  the  y 

n 

curve,  the  failure  is  not  going  to  affect  the  drift  mobility  calculation 
slgnlfloantly.  The  region  of  failure  is  too  small  to  have  any  signifi¬ 
cant  effect. 

The  analytical  formulas  whloh  are  derived  from  the  Perturbation 
Method  are  simpler  because  each  y*  can  be  calculated  separately'.  The 
aocuracy  of  each  yR  calculated  is  independent  of  the  number  of  equations 


in  (2.10).  The  formulas  for  tha  different  y“  given  in  (3,l3)-(3.20)  are 

n 

relatively  simple  and  easy  to  simplified.  Another  advantage  of  this 
method  is  that  only  the  first  two  terms  in  eaoh  expansion  are  needed  to 
provide  a  good  approximation  for  eaoh  yr.  The  resulting  analytical 
formulas  are  also  more  eompaot  and  simple. 


Onoe  y($,e  )  has  been  found,  the  eleotron  mobility  can  be  found  by 


the  relation  [141 


>> 


-  4h(2nkBT)1/2(e®-l) 

3icem<m*  )3/2(l/*.-  1/e,) 


(5.1) 


For  comparisons  between  the  various  method,  a  function,  G(o  ),  which  is 
defined  in  the  form  shown  in  (5.2)  will  be  used. 

am 

G (a  )•  f  y(?.e  )?  d$.  (5.2) 

This  funotlon  has  been  used  by  Fortini  et  al.  in  their  paper.  A 
comparison  between  the  factor  G(«  )  obtained  by  the  different  methods  is 
shown  in  figure  5.  The  values  of  some  of  the  G(<j  )  are  tabulated  in 
table  8.  Some  of  the  values  obtained  by  Fortini  et  al.  has  been  in¬ 
cluded  for  comparison.  The  number  of  equations  used  in  the  Iteration 
Method  for  obtaining  the  factor  G(e)  is  twelve.  In  the  Perturbation 

Method,  the  number  of  y™  summed  is  different  for  different  y  .  For  o»l, 

n  n 

eight  (i.e.  n-0  to  7)  terms  are  used  and  the  number  is  decreased  as  a 
lnoreases,  such  that  for  a  given  o  the  yR  used  are  yQ,  y^,  y^,... 

where  n-7-e  for  each  integer  a  *  The  number  of  terms  summed  for  each 
Yn  also  vary  with  n  as  described  in  section  4. 

From  table  &,  it  can  be  seen  that  the  results  from  the  Iteration 
Method  used  in  this  report  are  smaller  than  that  obtained  by  Fortini  et 
al  (141.  The  only  exception  is  when  <r  is  equal  to  one.  The  results 
obtained  in  this  report  are  smaller  beoause  fewer  equations  are  used 
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1.5413 

1.5246 

1.6173 

1.6190 

1.6186 

1.6784 

1.6584 

1.7322 


1.7301 


1.7953 


1.7745 


here  for  the  oaloulations  of  Q(e  ).  The  error  that  results  from  using 
fewer  equations  is  more  serious  for  smaller  e  .  The  results  Improve  as 
a  increases  and  the  error  is  less  than  one  percent  beyond  a  -2.5. 

This  is  also  an  indication  of  the  fewer  equations  needed  for  the 
Iteration  Method  when  a  is  large.  The  reason  why  the  result  obtained 
by  the  Iteration  Method  here  is  larger  than  the  result  of  Fortlni  et  al 
when  o-l,  oan  also  be  explained  in  the  same  manner.  If  one  looks 
at  the  plot  of  G(«  )  in  the  paper  of  Fortini  et  al  which  extended  the 
G (a  )  to  or  smaller  than  one,  one  notices  that  at~cr  -1,  G(o  )  is  near  a 
minimum.  Therefore  it  is  conceivable  that  when  not  enough  equations  are 
used  in  calculating  G (cr  ),  the  minimum  will  be  higher.  In  this  report  «r 
is  not  extended  to  values  smaller  than  one  because  the  mathematical 
approximation  is  poor  for  <r<l. 

The  Perturbation  Method  has  results  that  match  closely  those  from 
the  Iteration  Method  for  the  whole  range  of  a  except  for  the  small 
region  near  a  -1.  This  is  not  surprising  since  it  has  already  been 
shown  that  the  Perturbation  Method  works  best  for  larger  a  .  For  9 
greater  than  two  the  Perturbation  Method  gives  results  that  differ  from 
the  Iteration  Method  by  less  than  five  percent.  The  differences  drop  to 
less  than  1  percent  beyond  a  -3. 

The  first  Analytical  Method  in  (Seotlon  4.1)  is  worst  in  approxima¬ 
ting  the  factor  G(o  ).  At  a  -1  the  value  of  G(o  )  obtained  by  this 
method  is  about  half  of  the  expected  value.  However,  the  approximation 
becomes  better  as  a  increases,  and  the  best  approximation  is  near  o  -3. 
For  higher  a  ,  this  Analytical  Method  overestimates  the  funotlon  G(c  ) 
by  about  five  percent.  The  values  of  G(<r  )  have  the  largest  deviation 


from  the  expected  result  when  a  Is  snail  because  in  that  region,  the 
analytical  expressions  as  given  in  (4.20),  (4.22)  and  (4.23),  carry  too 
few  terms  to  make  a  good  approximation.  Moreover  for  smaller  a  ,  the 
convergence  of  the  yQ  expansions  is  slower  and  more  than  just  the  first 
three  terms  are  needed  in  each  expansion  for  a  better  approximation. 

The  analytical  formula  that  is  derived  from  the  Perturbation  Method 
(Section  4.2)  gives  a  good  approximation  for  the  whole  region  of  <x  . 

The  best  results  are  between  a  -2  and  3.5.  For  most  of  the  range  the 
analytical  formulas  over  estimate  the  factor  0(o  ).  However  on  the 
whole  the  result  matches  that  of  the  Iteration  Method's  results  rather 
well.  The  error  is  less  than  5  peroent  for  the  whole  range  shown.  The 
accuracy  that  is  found  for  the  region  where  a  is  small  is  unexpected. 
Since  the  analytical  formula  is  derived  from  the  perturbation  formula 
with  some  approximations,  it  is  not  expected  to  give  results  that  are 
more  aoourate  than  the  Perturbation  Method.  In  this  case  it  appears  that 
the  errors  from  the  Perturbation  Method  and  the  errors  from  the  approxi¬ 
mations  made  when  deriving  the  analytical  formula  compensate  each  other 
in  the  region  where  <r  la  small. 

After  the  value  of  G(e)  for  a  given  a  has  been  calculated,  the 
drift  mobility,  p,  at  that  temperature  can  be  calculated  by  substituting 
into  (5.1) .  The  drift  mobility  in  (5.1)  can  be  written  as 

U  -  -C  T1/2(ec  -l)G(e).  (5.3) 

m 

The  faotor  C  is  a  constant  whioh  depends  on  the  measured  parameter  m*. 

ZB 

•Mand  For  GaAs  Fortini  et  al  use  a  value  of  233  for  Cmbut  it  is 
found  that  with  the  values  of  the  parameters  given  by  them  (a*«0.073mo. 
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8*11.10  and  t-13.13),  the  value  for  C  is  227.  In  figure  S.  the  drift 

^  m 

nobility  for  GaAs  is  plotted  agalnat  the  temperature  with  0^*227.  As  in 
the  ease  of  the  plot  of  factor  G(cr  )  versus  a  .  the  analytical  formula 
derived  froa  the  Iteration  Method  shows  the  largest  deviation  froa  the 
Iteration  Method. 
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Figure  6.  The  drift  mobility  for  GaAs  calculated  from  the  four  methods 
are  plotted  against  the  temperature.  The  result  of  the  Iteration  Method 
is  denoted  by  the  solid  line,  the  result  of  the  Perturbation  Method  is 
denoted  by  the  dashed  line,  the  result  of  the  analytical  formulas 
derived  from  the  Iteration  Method  and  Perturbation  Method  are  shown  by 
the  dotted  line  and  dash-dotted  line  respectively. 
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6 .  ELECTRON  DISTRIBUTION  FUNCTIONS  FOR  HIGH  ELECTRIC  FIELDS 

For  high  electric  fields,  the  carrier  distribution  function  is 
found  by  solving  the  Boltzmann  equation.  The  validity  of  using  the 
Boltzmann  equation  for  high  electric  fields  has  not  been  fully  investi¬ 
gated.  However  ,  it  seems  for  the  range  of  the  fields  that  is  discussed 
in  this  report,  the  Boltzmann  equation  can  be  used. [25] 

In  III-V  compounds,  at  high  electric  fields  near  room  temperature, 
the  scattering  is  still  predominantly  polar  optical  in  nature  as  in  the 
low  field  case.  Other  forms  of  scattering  such  as  impurity  scattering  , 
space  oharge  scattering  and  acoustic  scattering  are  unimportant  and  will 
be  negleated.  For  electron  energies  that  are  greater  than  the  valley 
separation  between  the  central  and  next  higher  valley  in  the  conduction 
band,  intervalley  scattering  has  to  be  included.  In  this  report  inter¬ 
valley  scattering  is  assumed  to  be  unimportant.  Conwell  [26],  Conwell 
and  Vassell  [6]  have  discussed  the  various  approaches  that  can  be  taken 
when  intervalley  and  other  forms  of  scatterings  are  to  be  included. 

For  a  given  electric  field  $,  the  solution  of  the  Boltzmann 
equation  may  be  expanded  in  a  series  of  Legendre  polynomials,  Pn(cosd), 
where  &  la  the  angle  between  S  and  the  wave  vector  k.  The  distribution 
function  is  nearly  spherically  symmetric  if  the  fraction  of  the  carrier 
energy  lost  in  a  collision  is  small  and  if  the  scattering  is  not 
predominantly  in  a  given  direction  [27],  It  has  been  shown  that  under 
these  conditions  the  first  term  in  the  series  of  Legendre  polynomials  is 
dominant,  while  the  third  and  higher  order  terras  in  the  series  are  very 


small  [28],  Therefore  for  a  simple  model,  the  distribution  funotion  oan 
be  approximated  quite  accurately  by 

f(k)  -  tQ(t)  +  fjoosfr-  f0(£>  +  kgg(E)  (6.1) 

where  kg  is  the  component  of  ic  in  the  direction  of  fi. 

In  general  polar  optical  scattering  is  non-elastic.  For  III-V 
compounds  such  as  GaAs,  at  high  electric  fields,  the  carrier  energy  is 
much  larger  than  the  optical  phonon  energy  and  the  scattering  is 
reasonably  elastic  and  Isotropic.  In  most  cases  other  forms  of 
scattering  such  as  impurity,  acoustic  and  space  charge  scattering  are 
also  present.  Some  of  these  scattering  processes  are  elastic  and 
randomize  the  electrons  with  lower  energies  so  that  (6.1)  still  holds. 

In  this  report,  the  electron  distribution  funotion  is  assumed  to  to  be 
of  the  form  shown  in  (6.1).  Unlike  the  low  field  case,  the  spherically 
symmetric  part  of  (6.1)  is  not  the  equilibrium  Maxwell-Boltzmann 
distribution  function  but  is  a  function  in  which  the  heating  of  the 
electrons  by  the  strong  electrio  field  produces  an  elevated  effective 
electron  temperature. 

6.1  SOLUTION  FOR  HIGH  ELECTRIC  FIELDS  OF  THE  ELECTRON 

DISTRIBUTION  FUNCTIONS 

In  this  seotion  the  electron  distribution  function  in  high  electrio 
fields  for  polar  optioal  scattering  is  calculated  from  the  Boltzmann 
equation,  (1.1).  The  method  used  by  Conwell  and  Vassell  [6]  and 
Stratton  [5]  is  discussed  briefly.  In  this  method  the  difference- 
differential  equations  are  reduced  to  differential  equations  by  Taylor 


expansions  of  the  distribution  functions.  A  new  method  is  presented  in 
vhioh  the  ooupled  difference-differential  equations  for  the  symmetric 
and  asymsetric  parts  of  the  distribution  function  are  derived  by  expan¬ 
sions  in  1/E.  The  solution  from  both  methods  is  compared. 

The  collision  term  and  eleotrio  field  term  of  the  Boltzmann 
equation  for  an  arbitrary  energy  band  have  been  derived  by  Conwell  and 
Vassell  [61.  The  energy  band  is  assumed  to  be  parabollo  in  this 
section.  The  energy  of  the  carrier  electrons  is  assumed  to  be  much 
larger  than  the  longitudinal  optical  phonon  energy.  Since  (1.11), 

(1.12)  and  (1.16)  are  for  an  arbitrary  energy  band,  they  can  be  simpli¬ 
fied  and  substituted  into  the  Boltzmann  equation  for  the  case  of  a 
parabolic  energy  band.  For  parabolic  energy  bands  and  with  p-wave 
scattering  neglected,  0  is  equal  t,  and  o  and  o'are  zero.  Therefore 
(1.10),  (1.11)  and  (1.13)  reduoe  to  (6.2),  (6.3)  and  (6.4)  respectively: 


where  H(x-l)  is  the  unit-step  funatlon,  and  ^(x)  and  x  are  defined  as 
^(x)  -  2sinh"1(x)1/2 


(6.5) 


X  »  ?/o 


(6.6) 


and  fg (x)  does  not  have  the  sane  argument  as  fg  (£)  but  is  used  for 
simplicity.  Equations  (6.2)-(6.4)  have  been  derived  previously  (61. 


An  analytical  solution  for  the  eleotron  distribution  function  from 
equations  (6.2)— (6.4)  has  been  derived  by  Stratton  (51  and  Conwell  and 
Vassell  (61  for  x  >>  1.  They  solved  the  equations  by  expanding  fg(x+l) 
and  g(x+l)  in  Taylor  series.  This  transforms  the  ooupled  difference** 
differetlal  equations  into  differential  equations.  The  Inverse  hyperbo¬ 
lic  function  ,  ip(x),  is  expanded  in  terms  of  ln(4x)  and  a  series  in  1/x 
By  retaining  terms  up  to  second  order  for  the  Taylor  expansion  and  by 
keeping  terms  up  to  1/x2  in  the  expansion  for<jkx),  the  symmetrical 
part  and  the  asymmetrical  part  of  the  collision  term  are  reduced  to 
(5.61 

/dfjA  ,  -eEg  )d  ri(e‘  ♦l)ln(4x)df0]*f0ln(4x))  (6.7) 

\dt  ayn  /2o*h«x|dx  L2(e'1  -1)  dxj  ) 

jk^g  .  (6.8) 

These  equations  are  valid  only  for  very  large  x.  From  (6.8)  a  relaxa¬ 
tion  time  for  the  case  when  x>>1  can  be  written  and  is  given  by  (5,61 


.  -  «Bo  (*a  +1 
dt  fQ  a3y  V2m*h«x  \ea  -1 


-kgg  y(2m*h«x)  I  eq  -1  ^ 


(6.9) 


Equations  (6.4),  (6.7)  and  (6.8)  oan  be  substituted  into  the  steady 
state  Boltzmann  equation  to  solve  for  f  (x)  and  g(x).  Two  ooupled 


differential  equations  in  fQ(x)  and  g(x)  are  obtained  by  equating  the 
eoeffioients  of  the  first  two  Legendre  polynomials  to  zero  separately. 


-2eB  d(x3/2g) 
Tti/x  dx 


.  *Eo  Id 

/2a*hux( dx 


<* 

l(e +l)ln(4x)dfQ 
2(e®  -1)  dx 


+ 


fQln(4x)  j 


0 

(6.10) 


eJ^  l£o  +  -  0  (6.11) 

m*  dx  x 


Proa  (6.11)  g(x)  can  be  expressed  in  terms  of  dfQ/dx  and  substituted 
into  (6.10).  The  resulting  exact  second  order  differential  equation  can 
be  Integrated  onoe  to  give  a  first  order  differential  equation  in  fQ(x). 
This  first  order  differential  equation  is  found  to  be  separable  and  can 
be  integrated  to  give  [5,6] 


f 


0 


F(0)  exp 


ln(4x) 

+  dln(4x) 


(6.12) 


where  F(0)  is  a  constant  that  results  from  the  integration.  The  con¬ 
stants  b  and  d  are  defined  as 

b  -  4E2(eJ  -1)  (6.13) 

SEjfe*  +1) 

,  1  »ll  (6.14) 

d  "  2(^1) 

The  asymoetrical  part  g(x)  oan  be  found  by  differentiating  fQ(x)  and 
substituting  the  result  into  (6.11). 


There  are  two  main  draw  baoks  with  the  result  obtained.  The  most 
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important  is  that  fQ(x)  given  in  (6.12)  is  not  normalizable.  This 
problem  is  the  result  of  negleoting  intervalley  scattering.  Another 
drawback  of  (6.12)  is  that  it  oannot  be  extended  to  low  eleotron 
energies.  The  result  is  expeoted  to  be  valid  at  high  energies  because 
the  Taylor  expansions  that  are  made  for  fQ(x),  g(x),  and0(x)  are  for 
high  electron  energies.  Both  fQ(x)  and  g(x)  fail  around  x»0.25  when 
ln(4x)  is  zero.  For  x<0.25.  ln(4x)  becomes  negative.  In  fact  even 
when  x  is  slightly  larger  than  0.25.  a  point  is  reached  when  ln(4x)  is 
smaller  than  other  terms  in  the <£(x)  expansions  which  are  dropped. 
Therefore  calculated  values  of  tQ(x)  and  g(x)  have  failed  even  at  those 
points . 

In  this  section  an  attempt  is  made  to  extend  the  result  of  Stratton 

(51.  and  Conwell  and  7assellC6]  to  low  energies.  From  (6.12)  one 

notices  there  are  two  parameters.  6  and  EQ/E,  which  can  be  used  as 

expansion  parameters.  For  high  electron  energies,  expansions  of  b  and  d 

in  S  give  leading  terms  which  are  independent  of  5.  Therefore,  the 

expansion  in  &  is  not  expected  to  work  well  since  it  does  not  simplify 

(6.12)  much.  The  other  parameter.  Eg/E,  can  be  used  as  an  expansion 

parameter  since  it  is  smaller  than  unity  for  high  fields.  In  (6.11)  for 

very  large  E  the  term  with  coefficient  b  will  dominate.  The  coefficient 

b  can  be  brought  out  of  the  integral  and  the  exponent  can  be  expanded  in 

terms  of  1/b,  This  expansion  whioh  is  equivalent  to  an  expansion  in 
2 

(Eg/E)  appears  to  work  even  at  lower  energies.  Therefore  2E/3EQ  is 
chosen  as  an  expansion  parameter  and  is  represented  by  K. 

It  is  assumed  that  the  eleotron  distribution  functions  fQ(x)  and 
kg(x)  can  be  expanded  in  power  series  in  1/K: 
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f0(x)  "  £  *o(x)  («.1S) 

r.(x)  -  V(m)j  t{(x)  (6.i6) 

J 

where  f^(x)  is  used  to  represent  kg(x).  The  superscripts  on  fjj  and 
are  not  powers  but  are  there  to  indicate  the  order  of  the  terms  in  the 
expansions.  These  series  forms  of  fQ(x)  and  f^(x)  are  substituted  into 
(6.2)-(6.4).  The  resulting  equations  when  substituted  into  the  steady 
state  Boltzmann  equation  will  allow  us  to  find  two  coupled  equations  in 
fQ(x)  and  fj(x).  By  equating  the  coefficients  of  the  first  and  the 
second  Legendre  polynomials  to  zero  we  find  the  two  equations  are 
equation  (A4.8)  and  (A4.9)  in  Appendix  4.  From  them  equations  relating 
the  different  terms  in  the  expansions  of  fQ(x)  and  f1(x)  can  be  found  by 
equating  the  ooeffioients  of  different  powers  of  1/K.  The  first  few  of 
these  equations  are  shown  below: 

x  -  0  (6.17) 

dx 


d[xf?(x)]  -  0  (6.18) 

dx 


dtxfj (x)] 
3x  1 


xd(f0  (x)] 
dx 


-  N HfJcxJ-e®  fj|(x+l>|0(x> 


-H(x-l)^<x-l>  [fjtx-l)-.1’  f°<x)|j  «.19) 


,20) 


where  (6.17)  and  (6.18)  are  found  from  equating  the  ooeffioients  of 
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(1/K)°,  and  (6.19)  and  (6.20)  are  found  from  equating  the  ooeffloienta 
of  (1/K) l.  The  equations  found  from  equating  higher  powers  of  1/K  are 
given  in  Appendix  4.  The  various  terms  in  the  fQ(x)  and  f^(x)  expan¬ 
sions  oan  be  found  by  solving  these  equations.  The  functions  f0(x+l), 
fj(x+l)  and  4>(x)  are  not  expanded  into  Taylor  series  for  solving  the 
problem  because  the  main  interest  here  is  to  look  at  lower  electron 
energies  region  (i.e.  small  x). 

The  first  equation  to  be  solved  Is  (6.17).  It  is  Integrated  and 
the  result  is  f|j(x)  equal  to  a  oonstant.  This  result  is  not  surprising 
since  it  oan  be  deduced  from  (6.12).  For  large  x  and  large  E  the  term 
bx  in  (6.12)  will  dominate  the  denominator  of  the  exponent.  In  the 
limit  as  x  and  E  become  very  large,  the  exponential  is  reduced  to  one. 
and  fQ(  x)  approaches  F(0).  It  is  therefore  possible  to  conolude  that 
the  leading  term  of  the  fQ(x)  expansion  is  F(0), 

f°(x)-F(0).  (6.21) 

The  next  equation  to  be  solved  is  (6.18).  It  is  integrated  and  by 
taking  the  limit  as  x  approaches  zero,  it  can  be  seen  that  the  integra¬ 
tion  oonstant  has  to  be  zero  or  else  f®(x)  is  infinite  at  x-0, 

fj(x)  -  0.  (6.22) 

This  result  and  the  result  for  f®(x)  are  substituted  into  (6.19)  and 
(6.20)  to  solve  for  f*(x)  and  fj(x)  respectively.  To  integrate  (6.19), 
all  the  f®(x)  are  replaoed  by  the  oonstant  F(0).  This  first  order 
differential  equation  oan  be  reduced  to  an  equation  with  integral* of 
inverse  hyperbolic  sine  funotions.  The  integrations  are  oarrled  out  in 
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Appendix  4  and  the  result  Is  shown  below, 

f\(x)  -  ?SOUl=slJ  ((2x+l)sinh“Vx-/xTx+I) 

x  -  _  _ 

-H(x-l)I(2x-l)sinh  v7x=D-/xTx=IT]}  .  ($.23) 

Sinoe  fj(x)  is  zero,  integration  of  ($.20)  results  in  f*(x)  being  a 

oonstant.  The  constant  can  be  determined  by  matching  the  expansion  of 

fQ(x)  in  1/K  with  ($.12).  It  has  been  shown  that  for  large  E  and  large 

x  the  leading  term  in  ($.12)  for  an  expansion  in  1/b  is  F(0).  The  next 

2  2 

order  term  is  of  order  (EQ/E)  or  1/K  .  Sinoe  there  is  no  term  of  order 
&0/E,  the  1/K  term  in  ($.1$)  is  expected  to  vanish,  which  implies  that 
the  oonstant  fj(x)  is  zero. 

Next  the  coupled  equations  found  from  equating  coefficients  of 
(1/K)  are  Integrated  (  see  (A4.14)  and  (A4.1S)  in  the  Appendix  4).  By 
using  the  result  Just  found  (i.e.  f*-0)  the  right  side  of  (A4.14) 
vanishes  and  the  resulting  differential  equation  can  be  integrated.  By 
taking  the  limit  as  x  approaches  zero  it  is  found  ,  as  in  the  case  of 

i 

fj(x),  that  fj(x)  is  zero  or  else  the  distribution  function  is  infinite 

2  1 

at  x»0.  To  solve  for  fQ(x)  in  (A4.13),  the  result  for  fj(x)  obtained  in 

($.23)  is  needed.  The  resulting  equation  obtained  from  substituting 

($.23)  into  (A4.15)  oannot  be  integrated  easily.  Therefore  fg(x)  is 

written  in  terms  of  an  integral  of  f*(x)  as  shown  in  (A4.23).  The 

2 

oonstant  of  integration  ,  Fg(0),  is  included  so  that  the  final  expansion 

of  f q (x )  will  match  the  result  of  ($.12)  for  large  x.  This  integration 

oonstant  can  be  determined  by  looking  at  the  condition  on  fQ(x)  as  x 

approaches  •.  For  ($.13)  to  match  (6.12)  as  x— *  •,  f*(x)  needs  to 

2 

vanish.  Therefore  the  constant  of  integration  ,  FQ(0),  is  an  integral 


froa  zero  to  Infinity  of  some  funotions  of  f^(x),  sinh  (x)  and  x,  as 

2 

shown  in  (A2.24).  The  fin^l  fora  of  fQ(x)  is  shown  below. 


^(xJ-^dx^f^xJsinh”1/?  -  e®  f  j(x+l)|Qx+l)  slnh“Vx  -lj 

-H(x-l)  2e®  f}(x)sinh"Vx=I  -fj(x-l)L  (2x-l)  sinh 
|  1  1  f  7x75=1) 

-1  fdx  J2fj (x) sinh*"1/?  -  e®  f* (x+1)  f(2x+l)  sinh"1^  -ll 
3/T(  1  1  |/x7x+D  J 


(x)3inh_1/x  -  e®  fj 


e®  fj(x+l)f(ax±__ 
1  [/xTx+D 


♦H(x-l)|2e®  f*(x)sinh_1,/x=l-f*(x-l)  l(2x-i)  sinh-1/?11!  — 1 1 1  >  . 

I  1  1  [~x(x-l>  1 

(6.24) 

So  far  two  non-zero  terms  in  the  fQ(x)  expansion  have  been  obtained 
by  integrating  (A4.16).  The  result  of  the  integration  is  an  expression 
for  fj(x)  which  la  the  product  of  1/x  and  an  integral  with  f^(x)  as 
shown  in  (A4.26).  The  Integral  can  be  simplified  from  a  multiple  lntegr 
to  a  single  integral  by  using  integration  by  parts  as  shown  in 
Appendix  4.  The  resulting  expression  for  f*(x)  is  shown  in  (6.2S), 


fj(x)-l/x  j[f 


fj(x)-e®  fj(x+l)  1 1  (2x+l)sinh"1/x-/xT5t+D  1 


-H(x-l)  [f J(x-l)-e®  fj(x)  1 1  <2x+l)sinlT1/x=I-./x(x+i)] 
-  I  dxtClx+Dsinh^/x-yxfx+l)]  rdfg(x)-e®  df^x+l)! 

„  Idx  dx  J 

♦  /dxl(2x-l)sinh"1/x=l-/xTx=D]  dfj(x-l)-e®dfj(x)  ' 
J  dxu  dx°  j 

I  § 


I' 

(6.25) 


Even  though  f*(x)  written  in  this  fora  looks  more  complicated  than  in 
the  fora  it  had  before  the  integration  by  parts  was  done  (as  shown  in 
(A4.26) ) .  it  is  aotually  in  a  better  form  for  numerical  integration. 
Written  in  this  fora,  the  first  two  terms  can  be  calculated  by  straight 
substitutions  onoe  fJot>  has  been  oaloulated  froa  (6.24).  The  factors 


with  the  arguments,  x,x+l  and  x-1,  in  the  last  two  terms  are 
functions  of  f*(x)  as  shown  in  (A4.15).  Therefore  the  process  of  calcu¬ 
lating  f3(x)  and  f q (x)  involves  a  single  Integration  of  a  function  with 
f*(x).  In  this  report  it  is  assumed  that  the  first  two  terms  in  the 
expansion  (6.15)  and  (6.16)  are  good  enough  approximations  for  fQ(x)  and 
fj(x). 

6.2  ASYMMETRIC  PART  OP  THE  DISTRIBUTION  FUNCTIONS  FOR 
HIGH  ELECTRIC  FIELDS 

In  figures  7.  8  and  9,  the  asymmetrical  part  of  the  electron 
distribution  function.  f1(x)/F(0)  is  plotted  against  x  for  both  the 
methods  discussed  in  seotion  6.2.  For  the  method  used  by  Stratton[5] 
and  Conwell  and  Vassellt6],  fj  is  obtained  by  differentiating  (6.12)  and 
substituting  it  into  (6.11).  This  form  of  f^x)  obtained  from  (6.11)  is 
shown  in  (6.26). 

» 

f. (x)  -  -3K  x  dfft(x)  (6.26) 

1  (2N+1)  dx° 

ft(x)  -  fJ(N/K)  +  fJ(N/K)3.  (6.27) 

Equation  (6.27)  shows  the  formula  for  calculating  f^x)  that  is  derived 
from  the  new  expansion  method  suggested.  In  this  case  it  is  assumed 
that  the  first  two  non-zero  terms  in  the  1/K  expansion  of  f^(x)  are  good 


enough  to  approximate  it 


Figure  8.  Graphs  of  f,(x)/F(0)  from  the  two  different  methods  are 
plotted  against  x  forthe  ease  when  a  =1  and  K=5.  The  dashad  line  is 
for  the  method  of  Stratton/Conwell,  and  the  dotted  line  is  for  the  1/K 
expansion  method. 


For  figures  7  and  8,  the  value  of  a  used  Is  i  for  both  while  K  is  equal 
to  3  In  figure  7  and  K-5  in  figure  8.  For  figure  9  the  value  of  <r  and  K 
are  1.462  and  3  respectively.  The  value  of  1.462  is  used  for  a  in  this 
case  because  it  corresponds  to  the  oase  of  GaAs  at  room  temperature. 

The  fx(x)  obtained  by  the  1/K  expansion  method  is  valid  to  x-0. 

The  f^x)  obtained  from  differentiating  (6.14)  is  not  valid  around 
X-.2S.  There  is  a  very  noticeble  difference  between  the  result  obtained 
by  the  two  methods  for  x  less  than  2.  The  difference  becomes  less  as  x 
increases.  Beyond  x-9.5,  the  difference  is  less  than  5  peroent.  The 
largest  difference  ocour  around  the  maxima.  In  the  case  of  <r-l  and  K«3 
the  fj(x)  obtained  from  the  method  of  Stratton  [51  and  Conwell  and 
Vassell[6l  has  a  maximum  around  x».75  and  has  a  f1(x)/F(0)  value  of 
0.5607.  The  fj(x)/F(0)  obtained  by  the  perturbation  expansion  in  1/K 
has  its  maximum  around  x-1  and  has  a  value  of  0.4396.  Since  the  fQ(x) 
in  equation  (6.13)  is  found  by  assuming  that  x  is  very  much  larger  than 
the  optical  phonon  energy,  it  is  not  expeoted  to  work  well  for  small  x. 

The  fj(x)  that  is  derived  from  it  is  therefore  not  expected  to  be  very 
accurate  for  small  x  either.  Therefore  the  ourve  of  fj(x)/F(0)  which  is 
calculated  form  the  1/K  expansion  method  is  expected  to  give  a  better 
description  of  f^(x).  For  large  K  the  expansion  method  is  expected  to 
give  a  better  description  of  both  the  symmetrical  and  asymmetrical  part 
of  the  distribution  function  for  all  x.  This  is  true  because  the  series 
expansion  in  l/K  depends  on  K  and  not  on  x. 

The  spherioally  symmetric  parts  that  are  found  from  both  methods 
are  not  normalizable.  As  x  approaches  infinity  the  values  of  fQ  from 
both  the  methods  approaoh  the  oonstant  F(0)  and  do  not  go  to  zero.  The 
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reason  for  this  Is  that  for  large  energy  (large  x)  intervalley 
scattering  has  to  be  included.  In  order  to  account  for  intervalley 
scattering,  an  additional  term  has  to  be  added  to  the  collision  tera. 
With  intervalley  scattering  inoluded  in  the  oaloulations,  a  normalizable 
electron  distribution  funotion  is  expected. 

A  better  description  of  the  distribution  function  is  also  found  if 
a  non-parabolio  energy  band  is  included.  Conwell  and  Vassell[6]  have 
shown  that  the  effect  from  non-parabolloity  of  the  energy  band  is  not 
important  for  very  strong  electric  fields.  However,  the  inclusion  of 
non-parabolic lty  in  the  energy  band  is  of  interest  because  it  has  been 
claimed  that  in  some  compounds  with  large  separation  between  the  central 
valley  and  the  next  satellite  valley  of  the  energy  band,  negative 
differential  resistance  can  be  caused  by  polar  optical  scattering  in  the 
central  valley  alone  [22] . 


j*  •  *  -  «  -  •  ,  -  *  *  ,  “  .  » 
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7.  EFFECTS  OF  ENERGY  BAND  NON-PARABOLICITY  AND  P-WAVE  SCATTERING 


In  this  section  the  effeots  of  non-parabolioity  in  the  central 
valley  and  p-vave  scattering  on  the  electron  distribution  function  and 
the  drift  mobility  will  be  investigated.  Both  the  effects  have  been 
investigated  by  Hauser.  Glisson  and  Littlejohn [22] .  The  aim  of  this 
section  is  to  check  their  results.  Two  of  their  findings  are:  (i)  in 
certain  ternary  and  quarternary  III-V  semiconductor  compounds  the  nega¬ 
tive  differential  resistance  can  be  caused  by  the  non-parabolicity  in 
the  central  valley  alone,  (ii)  the  effects  of  p-wave  scattering  on  the 
drift  mobility  for  the  III-V  compounds  investigated  is  negligible. 

When  looking  at  the  effects  of  the  non-parabolioity  of  the  central 
valley.  Hauser  et  al(22]  assumed  an  t- k  relation  of  the  form  shown  in 
the  equation  below. 

(h*kl)/(2m*>  -  1(1  +  at)  (7.1) 

where  a  is  the  first  order  non-parabolic ity  factor  quoted  from  Conwell 
and  Vassell[6].  This  oonstant  is  different  for  different  compounds  and 
is  related  to  the  effective-mass  energy  gap  and  fc.p  interaction  matrix  P 
[9].  For  more  details  on  a  one  is  refered  to  reference  [6].  By  using 
the  energy  and  momentum  equations,  Hauser  et  al  [22]  were  able  to  show 
that  negative  resistance  can  occur  for  electron  scattering  in  the 
central  valley  alone.  The  drift  mobility  of  Al  -In  ,.As  was  calcu- 
lated  by  using  this  simple  model.  They  showed  that  the  results. from 
their  analytical  model  agree  very  well  with  another  set  of  results  that 
they  obtained  from  Monte  Carlo  calculations.  In  the  Monte  Carlo  calcu- 


la t ions  the  authors  made  two  sets  of  oaloulatlons:  one  that  Inoludes  p- 
wave  scattering  and  another  without  the  p-wave  scattering.  They  were 
able  to  show  that  the  effect  of  p-wave  scattering  is  neglgible. 

There  appears  to  be  a  oontradiotion  between  the  results  of  their 
simple  analytical  model  and  the  Monte  Carlo  oaloulatlons  in  relation  to 
the  effect  of  p-wave  scattering.  It  is  true  that  if  the  energy  band  is 
nearly  parabolio  then  p-wave  scattering  is  not  important.  However  for  a 
non-parabolic  energy  band,  the  p-wave  mixing  cannot  be  neglected  beoause 
non-parabolic  effeot  cannot  be  large  if  p-wave  mixing  is  small  [23]. 
Therefore  it  is  inconsistent  to  assume  a  non-parabolic  energy  band,  and 
no  p-wave  and  s-wave  mixing. 

It  is  the  aim  of  this  section  of  the  thesis  to  investigate  how 
significant  the  p-wave  mixing  is  in  a  non-parabolio  energy  band  of  the 
type  given  in  (7.1).  Calculations  are  aade  to  see  the  magnitude  of  the 
p-wave  contribution  to  the  asymmetrical  part  of  the  eleotron  distribu¬ 
tion  functions.  As  in  the  oase  of  parabolio  energy  bands,  the  steady 
state  Boltzmann  equation  will  be  used.  The  collision  term  in  (1.1)  has 
been  derived  in  Appendix  1.  and  is  given  by  (Al.10)  and  (Al.ll).  It  is 
assumed  that  the  scattering  due  to  the  s-type  wave  state  is  dominant. 

Therefdre  the  magnitude  of  o  and  o'  are  small  compared  to  a  and  a*  in 

2  2 

those  equations.  It  is  possible  to  drop  (oo')  in  A*  without  muah 

error.  A'  has  been  defined  by  equation  (1,12)  in  this  thesis.  Since 

al«l-o^,  all  the  terms  with  aa’  can  be  replaced  by  (l-o2) (l-o'2) .  This 

2  2 

results  in  some  terms  with  o  o9  whloh  oan  be  dropped »  The  symmetrical 

2  2 

part  of  the  distribution  function  after  terms  with  o  o'  have  been 
dropped  is  shown  in  equation  (7.2), 
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where  the  prime  in  1'  and  k*  impliea  that  the  faetora  involved  are  to  be 
evaluated  at  E+hu  aa  indicated.  The  aayometrical  part  of  the  distribu¬ 
tion  function  ia  given  by  equation  (AS. 3)  in  Appendix  S.  In  the  case 
for  whioh  there  ia  no  p-wave  scattering,  c-0  and  the  equation  reduces  to 
(A4.1)  for  a  non-parabolic  band.  It  is  assumed  that  the  symmetrical 
part  of  the  collision  term  can  be  written  as  the  sum  of  a  dominant  part 
associated  with  s-wave  scattering  and  a  perturbed  part  whioh  is  due  to 
p-wave  scattering 
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where  e  ia  a  smallness  parameter  which  will  be  set  to  unity  later.  In 

this  form  ,  the  s-wave  associated  dominant  part  represents  those  terms 

2  2 
without  o  ,  while  (df_/dt)  „„  consists  of  those  terms  with  o  . 
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When  using  the  collision  term  in  the  Boltzmann  equation,  it  is 
convenient  to  express  the  k  and  k’  in  terms  of  the  electron  energy  . 
For  the  form  of  the  non-parabolio  energy  band  in  (7.1).  approximations 
are  made  in  order  to  arrive  at  a  term  whioh  is  considered  simple  enough 
to  be  used  later  on.  The  assumption  used  in  making  the  approximations 
is  &>>he.  This  assumption  is  appropriate  beoause  a  high  eleotrlo 


field,  E,  leads  to  a  high  average  eleotron  energy.  Moreover  it  appears 
that  p-wave  scattering  is  more  pronounoed  at  higher  electron  energies. 
The  three  approx iaat ions  that  are  made  involve  the  simplification  of  the 
term  (lc*  +k  )/kk',  the  expansion  of  fQ(t)  into  Taylor  series  and  the 
dropping  the  lnl(k*+k)/(k'-k)1  terms.  The  details  of  the  approximations 
are  in  Appendix  5.  The  simplified  form  of  the  symmetrical  part  of  the 
collision  term  due  to  p-wave  scattering  is  shown  in  (7.4). 


/dfjA  _-2eEno2  j4dh«f0<£)  +  hMdfgtOTl+la  +(2N+l)2ah»'| 

\dt  I  '  We  d t  ^ 


'c  sym 
p-wave 


+(h<*)2d2fft(t)  f"(l+2o£)  (2N+l)+2ahuil 

2  df  J 


(7.4) 


Prom  the  steady  state  Boltzmann  equation,  (1.1),  the  rate  of 
ohange  of  the  symmetrical  part  of  th  eleotron  distribution  function  in 
the  lower  valley  oan  be  written  in  the  form  shown  in  (7.5)  by  assuming 
no  intervalley  sacttering 


2e2E2 

3m*(l+2at) 
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[ JS* P Lff^T  _L£o  1  +  /!£o\  +  c/d_fo\  =0 
16+et  J  d£ll  hw  J  l+2at  J  \dt  /fl_9u\dt  / 


(7.5) 

where  the  non-parabolloity  in  the  energy  band  of  the  form  in  (7.1)  has 
been  used  for  the  field  term.  The  collision  term  has  also  been  written 
as  the  sum  of  a  s-wave  scattering  and  a  p-wave  scattering  terms  as  given 
in  (7.3).  The  only  difference  between  (7.5)  and  the  equation  in 
reference  [6]  is  the  extra  term  that  has  been  added  to  acoount  for  the 
p-wave  scattering,  t  is  the  polar  optioal  relaxation  time  for  £>>h<«  as 
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dlsoussed  in  section  6.  For  the  non-parabollolty  considered  in  this 


section,  it  is  given  by  [61 


11/2  e*  -  1  (t*al2)1/2 


r2a^h«1  e  -  1 
l  eE*  I  e®  +  1 


(7.6) 


r(hw)  (l+2at) 


The  collision  term  where  no  p-weve  scattering  is  assumed  has  been 
derived  for  a  non-parabolic  band  by  Stenflo  [29]  and  Conwell  et  al  [61, 
Por  convenience  the  form  derived  by  Stenflo  [291  will  be  used  here. 
With  the  expression  from  (7.4)  and  (7.6),  (7.5)  oan  be  written  in  the 
form  shown  below 
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(l+2at)(2N+l)+2dh« 


(7.7) 


It  is  assumed  that  the  solution  of  (7.7)  oan  be  written  as  the  sum  of  a 
dominant  term  f®(x)  due  to  s-wave  scattering  and  a  small  perturbed  term 
ef2(x)  caused  by  p-wave  scattering. 
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(7.8) 


where  i  is  a  smallness  parameter  to  be  set  to  one  later.  When  this  sum 
is  substituted  into  (7.7),  the  terms  in  the  equation  oan  be  reordered  in 
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teras  of  «°  and  «2.  The  eoefflolenta  of  each  power  of  s  oan  be  set  to 
zero  and  two  equations  are  obtained: 


2e2B2 
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(l+2a£)d  £/lft(d  d£  I  lh<i»(l+2a£  )J) 


Wc(l+2a£)d 


_2eE(,O2r4gft(jf0  +  ft<#dfJ(1+2ae>  +  (tjjjj)2  l!fo(l+2<£)(2M+l)l  A 

-hie  l  0  df  2  dtz  r  °* 

(7.10) 

The  first  of  these  equations  is  the  one  that  is  obtained  by  equating  the 
ooeffloients  of  the  first  Legendre  polynomial  PQ  to  zero  in  the 
Boltzmann  equation  with  p-wave  scattering  neglected.  This  equation  has 
been  solved  by  Conwell  et  al  [6]  for  an  arbitrary  energy  band.  The 
solution  that  they  derived  is  shown  below. 


f“  -  F(0)  exp-J 


pJ  (V)4  ln(4n/q»>  dx 
/  lb2n  +  d(n')4in(4n/n’) 


(7.11) 


where  xmUhu  is  defined  in  seotion  6,  b  and  d  are  defined  by  equations 
(5.14)  and  (5.15)  in  section  5  and  q  is  the  energy  function  0  normalised 
by  the  hw.  For  the  type  of  non-parabolioity'oonsidered  here. 


q-(£+a£.  )/(7ta»)  and  the  prime  In  q  indicates  differentiation  with  respect 
to  x.  F(0)  is  a  oonstant  of  integration.  In  Conwell  and  Vassell's  [6] 
paper,  F(0)  has  been  omitted.  For  a  non-parabolic  energy  band  (i.e. 
a«0),  (7.11)  would  reduoe  to  (6.13)  which  is  the  symmetrical  part  of  the 
electron  distribution  function  when  there  is  no  mixing  of  p-  and  s-wave 
functions. 


To  see  how  much  p-wave  scattering  affects  the  electron  distribution 

function,  equation  (7.10)  has  to  be  solved  for  f*(x).  There  are  three 

2 

terms  within  the  bracket  that  are  multiplied  by  the  factor  o  .  The  last 

of  these  three  terms  is  a  produo t  of  the  second  derivative  of  f|J(x)  and 
2 

(hu)  ,  which  1s  small  compared  to  the  other  two  terms  and  is  dropped. 

It  has  also  been  found  to  be  convenient  to  work  with  the  dimensionless 
variables  x  and  q  whioh  are  introduced  in  (7.11).  Equation  (7.10)  can 
be  Integrated  once  to  give  a  first  order  differential  equation  which 
still  contains  an  Integral  term  as  shown  below 

-0  ♦  „  n'* _ 

dx  bq  -Kin*  ln(4q/q')  bq  +dq'  ln(4q/q' 

2 

In  this  equation,  c  is  also  a  function  of  £,  and  is  written  in  terms  of 
x  so  that  the  right  hand  side  can  be  Integrated.  From  (1.4)  and  (1.5), 
it  is  possible  to  obtain  an  expression  for  o  in  terns  of  and 

Only  Eq  is  a  function  of  k  and  is  given  by  £Q«<h  k  /(2mQ),  where  mQ  is 
the  mass  of  electron  in  vacuum.  Since  m^  and  the  effective  mass,m*  are 
related  through  mQ»m*(l+£p/£g),  c  oan  be  expressed  in  terms  of 
1)V/(2m*),  £p,  and  £g.  This  expression  can  then  be  changed  to  into  one 

2  2  7. 

that  is  written  in  terms  of  q’  through  the  relations  71  k  /(2m •)-£+«£ 


/2c2d[q,4fJ  ]  dx  . 

*  (7.12) 


and  a  -£^/[£g(fg+^)*] .  The  various  substitutions  mentioned  above  have 

2 

been  oarried  out  In  Appendix  5,  and  the  final  result  for  o  is  a  simple 
function  of  V. 

o2-|(l-l/i,»)  (7.13) 

This  expression  is  substituted  into  the  right  hand  side  of  (7.12)  and 
integration  by  parts  is  used.  The  right  hand  side  of  (7.12)  can  also  be 
simplified  a  little  by  expressing  the  coefficients  of  fj  in  terms  of  the 
derivative  of  log  [F(0)/f®]  from  the  result  obtained  earlier  for  f®  in 
(7.11).  Therefore  (7.12)  oan  be  rewritten  as 


>r 

df  J+f Jd[in(F (0)  /f  J)  ] -  V2 _ U  '2f jj-f jjn  '-2ah«J  f°  dx] 

dx  dx  bn2+dn'*ln(4W)  * 

(7.14) 

It  is  possible  to  transform  the  right  hand  side  of  the  above  equation 
into  the  form  of  an  exact  integral  by  the  integrating  factor,  F(0)/fjj. 
With  the  integrating  factor,  (7.14)  can  be  integrated  from  a  lower 
limit,  X,  to  x.  This  lower  limit  la  ohosen  suoh  that  f2/f®  vanishes  at 
x-A.  The  lower  limit  is  found  to  be  one.  After  the  integration,  F(0) 
oanoels  out  from  both  side  of  the  equation,  and  the  result  of  fQ(x) 
expressed  as  a  ratio  of  fj  is 


The  faotor  F(0) (1-exp (1/d) )  is  an  approximation  for  the  integral  from 
aero  to  1  of  f^.  This  has  been  used  to  overoome  the  singularity  that  fQ 
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has  at  x«0.25,  as  diseussed  In  seotlon  6. 


The  effeot  of  the  p-wave  scattering  on  the  mobility  of  electrons 
can  be  seen  by  looking  at  the  asymmetrical  part  of  the  electron  distri¬ 


bution  function  g,  whioh  is  expressed  in  terms  of  fQ  from  (4.12). 
However  the  differenoe  here  is  that  in  this  case  fQ  oonsists  of  a  term 
associated  with  s-wave  scattering  and  an  extra  term  associated  with  p- 


wave  scattering. 


•^heEt 

a*hwi) 


•  dx  dx 


(7.16) 


where  the  first  term  is  due  to  s-wave  scattering  and  the  term  with  e  is 
due  to  p-wave  scattering.  P  is  the  expression  within  the  integral  in 


(7.15) 


b^  +da  *4ln(4n/n 


x 

_  n'2-n,-2ohw^P(0«l-e1/d)+  J  fjdxj 
.  i  . 


Since  the  mobility  is  proportional  to  the  integral  of  the  asymmetri¬ 
cal  part  of  the  electron  distribution  funotion,  the  effects  of  p-wave 
scattering  on  the  mobility  can  be  seen  by  looking  at  how  much  it  affects 


the  distribution  function. 


7.1  RESULTS  WITH  P-WAVE  SCATTERING  INCLUDED. 

The  ratio  fj/f®  as  given  in  (7.15)  has  been  numerically  integrated 
for  GaAs  and  A1  In  ^As.  The  results  are  plotted  in  figures  10  and  11 
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Figure  10.  Plots  of  the  normalised  symmetrical  part  of  the  distribution 
function  due  to  p— wave  scattering,  t^/FO  versus  x  for  GaAs.  Two 
different  values  for  E/E_,  3  and  6,  are  used  here  as  indicated.  The 


distribution  function,  f*  against 


The  integrations  have  been  carried  out  on  an  Apple  II  oomputer  with  a 
simple  numerical  integration  program.  For  eaoh  oompound  two  values  of 


E/Eg  are  used  in  order  to  see  the  variation  of  fg  with  eleotrio  field. 
The  values  of  a  have  been  calculated  from  the  known  values  of^  and 
for  eaoh  oompound.  The  range  of  the  dimensionless  energy,  x,  is  chosen 
suoh  that  it  is  about  equal  to  the  difference  in  the  energy  of  the 
central  valley  minimum  and  the  next  higher  valley  minimum.  Therefore, 
the  effect  of  p-wave  scattering  is  studied  for  the  energies  below  that 
for  whioh  the  transferred-eleotron  effect  becomes  Important. 


In  Figure  10  is  the  result  for  GaAs.  At  room  temperature  'hw  is 
about  0.03(eV  and  the  valley  separation  between  the  minima  of  the 
central  valley,  (D  and  that  of  the  next  higher  valley (l.)  is  about 
0.31eV.  Therefore  for  x>10.  Intervalley  scattering  has  to  be  Included 
in  the  collision  term.  For  a  given  electron  energy  (x),  the  effeot  of 
p-wave  scattering  decreases  as  the  applied  field,!,  is  increased.  For  a 
given  applied  field  the  effects  of  the  p-wave  scattering  increases  as 
the  energy  of  the  electron  increases.  The  increase  is  more  rapid  for 
lower  electric  fields.  However  the  perturbed  part  is  still  very  small 
for  the  whole  range  of  x  investigated  since  fg  is  always  less  than  0.3 
peroent  of  fg.  The  contribution  of  the  p-wave  scattering  to  the  asymme¬ 
trical  part  of  the  distribution  funotion  is  shown  in  figures  11  and  12. 
In  these  figures,  the  asymmetrical  part  f1  and  its  s-wave  component  and 
p-wave  component  are  plotted. 


As  can  be  seen,  the  p-wave  contribution  to  fj,  hence  the  mobility,  is 
small  for  GaAs.  Even  at  the  high  energy  end  where  x-10,  the  p-wave 
accounts  for  only  10  percent  of  the  scattered  electrons  at  that  energy. 
Therefore  the  contribution  of  p-wave  scattering  to  the  mobility  is 
expected  to  be  very  much  less  than  10  percent  since  there  are  more 
electrons  at  the  lower  energies  portion  of  the  distribution  function 
which  is  less  affected  by  the  p-wave  scattering. 

In  A1  2Jin  75A3,  the  effect  of  p-wave  is  more  important  as  can  be 
seen  in  figure  10.  The  size  of  f*  varies  from  0  to  3  percent  of  f®, 
which  is  larger  than  the  case  for  GaAs.  However  the  effect  is  still 
small.  As  in  the  case  of  GaAs,  f^  increases  with  increasing  x  and 
decreases  with  increasing  E.  The  fraction  of  the  asymmetrical  part  of 
the  distribution  function  due  to  p-wave  scattering  in  shown  in  figure  13 
and  14  for  E/EQ»3  and  6  respectively.  The  effect  of  the  p-wave 
scattering  is  more  important  here  than  in  the  case  of  GaAs.  At  x«30, 
about  a  third  of  the  electrons  at  that  energy  are  scattered  by  p-wave. 
However,  as  discussed  earlier,  the  effects  on  the  mobility  is  small 
since  there  are  more  electrons  in  the  lower  energy  portion  of  the 
distribution  function. 

7.2  DISCUSSION  OF  THE  RESULTS  WITH  P-WAVE  SCATTERING  INCLUDED 

Various  approximations  have  been  made  in  deriving  (7. IS)  and 
(7.10.  These  give  an  estimate  of  the  effect  on  the  electron 
distribution  function  due  to  p-wave  scattering.  The  conclusion  is  that 
p-wave  mixing  does  not  have  a  significant  effect  on  the  electron  distri¬ 
bution  funotion  and  hence  the  drift  mobility.  It  can  be  seen  from 
figures  11-14  that  the  contribution  to  the  asymmetrical  part  of  the 


electron  distribution  function  is  less  than  5  percent.  Since  'the  p- 
wave  scattering  does  not  affeot  the  drift  nobility  significantly,  the 
effeot  of  the  non-parabolloity  in  the  band  structure  oan  not  be  as 
laportant  as  has  been  claimed  by  Hauser  et  al  [22].  Their  suggestion 
that  the  negative  differential  mobility  in  oertain  ternary  and 
quarternary  III-V  compounds  occurs  for  scattering  within  the  central 
valley  only  if  the  energy  band  is  non-parabolio  band  may  be  wrong. 

Rather,  the  calculations  for  Al  25in>7SAs  seem  to  indicate  that  the 
negative  differential  resistance  for  elec trio  fields  near  the  peak  for 
the  velocity-electric  field  curve  arises  from  electron  scattering  in  the 
central  valley  of  a  parabolic  energy  band.  This  conclusion  is  reached 
not  only  based  on  the  findings  here  that  the  effect  of  p-wave  scattering 
on  the  eleotron  distribution  is  negligible  but  from  their  Monte  Carlo 
oaloulatlons  that  p-wave  scattering,  hardly  alters  the  velocity-electric 
field  ourve. 

To  understand  the  possibility  of  negative  differential  resistance 
in  III-V  semiconductor  compounds  due  to  scattering  in  the  central  valley 
alone  for  a  parabolio  energy  band  one  oan  oompare  the  results  for  QaAs 
and  Al  2^in  As.  Por  GaAs  the  peak  of  the  veloeity-eleotrio  field 
ourve  occurs  at  a  value  muoh  smaller  than  E-3BQ/2«9.3  kV/cm  whereas  for 
Al  23In  7JAs  the  peak  of  the  veloolty-eleotric  field  ourve  occurs  at 
roughly  E*3E0/2«5.2  kV/om.  The  value  of  E-3Eg/2  is  significant  since 
from  the  equation  for  polar  optioal  scattering  we  have  found  for  K<<1  an 
expansion  in  K  is  suitable  and  for  K>>1  an  expansion  in  1/K  is  suitable 
for  the  perturbed  distribution  function  that  determines  the  drift 
mobility.  From  the  fora  of  the  two  series  one  oan  reasonably  expect  a 


maximum  in  the  drift  mobility  for  K-l.  For  GaAs  the  energy  gap  between 
the  central  valley  and  the  closest  lying  upper  valley  of  the  conduction 
band  is  relatively  small.  Therefore  the  maximum  is  not  reached  at  E-9.3 
kV/om  because  a  significant  number  of  electrons  are  transfered  to  the 
upper  valley  at  much  lower  electric  fields.  On  the  other  hand,  for 
A1  2jin  As,  the  energy  separation  between  the  central  valley  and  the 
closest  upper  valley  is  considerably  larger  than  it  is  in  GaAs  and 
additionally  the  physical  parameters  of  A1  2jln  7JAs  give  a  smaller 
value  of  3Eq/2  (5.2  kV/cm) .  Hence,  for  A1  2JIn  7JAs  the  peak  in  the 
velocity-electric  field  curve  occurs  for  an  electric  field  such  that 
few  electrons  are  tranferred  to  the  upper  valley. 

For  the  physical  mechanism  that  causes  the  negative  differential 
resistance,  consider  the  hypothesis  that  the  electrons  are  heated  by  the 
electric  field  to  the  extent  that  the  increase  in  the  drift  mobility 
with  increasing  electric  field  is  offset  by  the  decrease  in  the  drift 
mobility  with  increasing  effective  electron  temperature.  From  (7.15) 
and  (7.1<S)  one  can  verify  that  it  seems  plausible  che  denominator  in 
the  integral  is  effectively  a  temperature  that  increases  as  the  electric 


field  increases 
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9.1  APPENDIX  It-  SYMMETRICAL  AND  ASYMMETRICAL  PARTS  OP  THE  ELECTRON 
DISTRIBUTION  FUNCTION. 


The  oolllslon  operator  for  a  non-parabolic  band  is  given  by 

k  +l(N-llH'!N)|2L(1t,!c')8(£Jj-^l-D«)f(ie') 

-I  (N+l  I H '  I N)  1 2L  He  ,le ' )  5  , -^+h« )  f  (S ) 

-I  (N-llH'|N)  |\<!c.lc*)S(2g,-££-D«)f  (£)]  (Al.l) 

The  sunmation  can  be  changed  to  an  integral  over  k-spaee  in  spherical 
coordinates.  The  scattering  matrix  is  given  in  section  1. 

1C  "  ygg*  j  |  J  Js^g-^.+R^tfdcMe^-frtc)] 

*  0  °+8(€^,-^-n*.)[f(C*)-e®f(ic)^.£|*£lj|»2sin  d*cifdk» 

(A1.2) 

Since  the  overlap  integral  is  given  by  LOc.lc'Waa ’+00*0031*) 2  where O' is 
the  angle  between  £  and  £',  the  rate  of  ohange  of  the  distribution 
function  oan  be  written  as 

1C  *  sjjjg!  fff  {*<V‘ic*4,ta,tf<k,,*''f<k,i  j 

’  ‘  If (k-)-.<,f(k)tf^2S^L!i:  sm«4*dfaic 

K  *  'fc»Z+lrZ-21 rVCna#' 


(A1.3) 


It  is  assumed  that  the  distribution  function  f(k)  can  be  written  in  the 
form  f(£)-f0(t)  +kgg(€)  where  kg  is  the  component  of  ic  in  the  direction 
of  E.  If '(Ms  the  angle  between  ie  and  I,  f(tc’)  can  then  be  written  as 


f(k,)*f0(a-,)+k,g{0*)  Csin  0  3in*eos(f-^  +  oos^eos^]. 


(A1.4) 


f  and  ^  are  the  azimuthal  angles  of  the  electric  field  and  k*  vector 
respectively.  The  integration  of f over  2n  results  in  the  term  with 
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ooslf-t)  equals  to  zero. 

If-  fteEnN  Jj  *C«#>-r0C«-lcB*C«I 


dt 


r  o 


+6(£Jc"^»  •"hw)  tf (t ' ) +k  'kEoos  g(2')-e<,f0(£)-efflcBg(e)n 

(..•««  WjV1  d(HM,rtdk. 

(A1.5) 


k'2  +k‘-2kk'eosd’ 


H  — 


|(k  '  +k ) 2-  (k  '  -k )  2-ln  I  k ' +k  I 


Integration  of  from  0  to  it  gives 

iSSjUl  {» (tj-tg ,♦»»>  [f„ ( V )e’-f0(t) -kggCt)] 

•  te,  H.ISSI-, «-^,-isam 

+(ce  • > 2  J  4-2  (JLlI+ld} 

(  2(2kk*)  ( 2kk  * ) 

+  l?(*lc*V  +tl“)  «^l'kgg(t»)/k  +  dCf^-^.+nwJk'kggtSM/k] 

«-(aa')2Ll  -  2 CkJL^ikllln |]£l+kj"}  -  (oc»)2 
|kk'  (2kk')z  Ik'-klJ  3kk' 

+(oc*)2(k,2+k2)  r-«c»k3-«rif»3  +2(k,2+k2)ln|k^+k/| 
2kk»  l 

'co '  J^gkkl 


k^klj) 


2(2kk') 


+2aa* 


l=ak!lC.+20c '2+k2>  lnlk '+k 


2(2kk'p 


k'-kl 
k,2dk' 


Ik'-k 


(Al.«) 


A  k*  in  the  denominator  oan  be  factored  and  oancels  part  of  k'  in  the 
numerator .  Then 


d£  -  heEnW  ^d(k'2)n8(6^,+he)  [fQ(  A'^-fgtO-kggCt)] 

\  +» (t^ , -he)  [f Q < V  >  (* )  -e^kgg  ( £)]} 

|A,2ln|k*4k|-A'oc*-aa'co'j 


+kg/k  [6  <t£tg ,  +ft»)  e*g  (t' )  +8  (t  ^ , --ft»)g  ( t)| 

... 2.  -l' 


T.  .2  .  2. 


.  .2.  . 


.1 


9.1  APPENDIX  1:-  SYMMETRICAL  AND  ASYMMETRICAL  PARTS  OP  THE  ELECTRON 
DISTRIBUTION  FUNCTION. 


The  collision  operator  for  a  non-parabolic  band  Is  given  by 


^2[l(N+llH»lN)l2L(lc,Ie')A(^-«Rt+he)f(Ic*) 
k  +1  (N-llH'iN)  l2L(lc,Ic,)8(£g-8|*,-'hw)f  (£') 

-I  (N+llH*  |N)  1 2Lflt.Tc * )8 ,-tg+he)  f  (Jt) 

-I  (N-llH'iN)  !2L(lt,le')&(Ejg,-g|.-N«)f  (lc)]  (Al.l) 

The  summation  can  be  changed  to  an  integral  over  k-space  in  spherical 
coordinates.  The  scattering  matrix  is  given  in  section  1. 


df 

dt 


tteE-N 

inrJ* 


fjj  Cf (C *  )e*-f  (lc)  1 

*  °  °+8(«t,-eft-!ta.)Cf(Ic*)-e<rf(k)^Lj|^ll^'2sin  define’ 


(A1.2) 


Since  the  overlap  integral  is  given  by  LOc.fc’J-taa'+oe'eosU)2  where is 
the  angle  between  ft  and  1c* .  the  rate  of  change  of  the  distribution 


function  oan  be  written  as 

g  -5^  fj 

♦8(L.-g.,-hw)tf(kM-e<Tf(k)ii(-lax~q.,Coa^  k*  sinOdOdf&k 
K  '  k’2+k  -2kk'Cos*; 

(A1.3) 


It  is  assumed  that  the  distribution  function  f()c)  can  be  written  in  the 
form  f(1c)«f0(E)  +kgg(t)  where  kg  is  the  component  of  ic  in  the  direction 
of  $.  If  <9is  the  angle  between  k  and  f(tc')  oan  then  be  written  as 


fCk'J-fQCe-M+k'gt^Olsin^ain^cosff-^  +  eos^oos^l.  (A1.4) 

f  and  'p'  are  the  azimuthal  angles  of  the  electric  field  and  k'  vector 
respectively.  The  integration  of  f over  2»  results  in  the  term  with 


where  A'-aa'+k'+lj.2  eo*.  It  is  assumed  that  theC-k  relation 
2kk’ 

is  given  by 


hV  -  n 

2m* 

with  d(k*2)  -  2g*  dn  d£'  .  we  find 
h*  d£' 


1C " 

*»«{«,;  ,-h.)  [f„  (  f  )  ,1)  -,\e  (*>]}• 

|A*2lnlOc*+k)/(k'-k)!-2A'cc’-aa'ec »J 

+k1!/k  [8 ea« (£)  +s (€  fa  ,-H»)g(t'5] . 


k,2+k2  A*ln|k^+k|-A,2k'-(oo ')2k' 


k  2  Ik'-k 


(Al.8) 


This  when  Integrated  gives 


df-®50 W  1 1 [fo( c * *<’Lf 0 (t)  "kE*  f «]  [* ' 2 In  I  (k '  +k )  /  (k '  -k )  I  -Aoc  '-aa'ceiJ 

'  +kgeVt)  /(k,2+k2)  A'ln|k^+k|-A’2k,-(ee')k'2\{  do 
\  k  2  (k'-kf  “3”  /(<£' 

+|[f0(4’)-e<,f0(d)-efflkBg(tj  [A'2ln!  (k»+k)/(k’-k)  l-A'eo’-aa'cc'] 

+kgg(t,)/k,2+k2  A^nlk’+kl  -A,2k’-(eel)V\  |d(l  ) 

Hr  \  k  2  !k'+kf  3  lidf  (  (A1.9) 


This  rate  of  change  of  the  distribution  funotion  due  to  collisions  oan 
be  separated  into  symmetrical  and  asymmetrical  parts  where  the  symmetri¬ 
cal  part  is  given  by 


I  /  (f  Vi 

2 

\ 

/  df ' 
\dt> 

1  -*B 
c-sym  ^ 

(f0(f>ew-f0(t))| 

i 

A'zln 

k'+k 

k'-k 

-A'eo’-aa'ee'j 

C 

r-£** 


+  |f0(tM^fff0(£)j|A'2ln|Stl+Si|  -A'oe’-aa'eo'j  Oft'  ) 


at'  / 

(A1.10) 


and  the  asyaanetrical  part  is 

/ df\  -  eEAMk0  (  e®g(£')/lc*2+lc2  A'lnlk'+kl-A'kMoe')2^ 
We-asy  *  j  -  \  “TT  T  IkM?!  T“  / 

-  g(«)|A'2lnl(k*+k)/(k'-k)|  -Arcc'-aa'ec'J  dfl» 


t'-fc+ft* 


+  g( &*)/k,2+k2  A;  lnjkj+kl-A'kMccoV) 
k  \  k  2  Ik'-kl  / 

-e®g (£>  (  A '  2lnl  (k  '+k)  /  (k  *-k)  I  -A 'cc '-aa 'cc  dO ' 

V  /  dt' 


1' 

(Al.ll) 


9.2.  APPENDIX  2:-  PERTURBATION  METHOD 


In  the  Perturbation  Method,  e_ff« 6  ia  used  as  the  expansion  parame¬ 
ter.  The  yQ  in  equation  (2.10)  is  expanded  in  the  form  shown  below. 

vj&  ->1 ,J  <“•» 

for  n>l,2,3... 

The  superscripts  of  yjj  indicates  the  order  of  the~terms  in  each  expan¬ 
sion.  The  coefficients  in  (2.10)  are  also  written  in  terms  of  S. 

Therefore  A  ,  B  .  C,  and  D  become  A®,  B^+B^/S,  C”1  and  Dn6n  respect! 
nnn  n  nnn  n  n 

vely.  Each  superscript  indicates  the  power  of  S  by  which  the  given 
coefficient  la  multiplied.  In  terms  of  the  new  representations,  (2.10) 
can  be  written  as 


J  «  j 


A1  E B®  £y^J  ♦  B^/8  £y^5J  ♦  C^/6  £y^ 


-dJs 


(A2.2) 


A°  +BU  ^n8^1  *  D“5* 

j  n"  j  °  6°  j  ”  “ 


By  equating  the  ooeffloients  of  the  different  powers  of  6,  the  set  of 
equations  obtained  from  the  coefficients  of  6-1  is  given  ln(A2.3X 


<T1y° 
C0  Y0 


•to 


—  a 

*  T. 


♦  C 


-1  0 


-  0 

-0 


B2  Y2  +  C2  Y2 


B-V‘  -0 

n  rn 


(A2.3) 


It  has  been  assumed  that  the  infinite  set  of  equations  in  (A2.2)  is 
trunoated  at  the  nth  equation.  The  result  from  solving  (A2.3)  is 


■  0  for  n>l. 
n 


(A2.4) 


The  next  set  of  equations  can  be  found  by  equating  the  coefficients  of 
5°  in  (A2.2) . 


*  c:lr}  «D? 
O'O  0  rl  0 


(A2.5) 


,0  0  _0  0-1  1  ‘ 

Vn-l+BnVBn  Yn 


When  (A2.4)  is  substituted  into  (A2.5),  the  set  of  equations  is  reduced 
to  Just  two  equations,  (A2.7a)  and  (A2.7b).  The  general  result  form 


solving  (A2.5)  is 


y  ■  0  for  n>2 
n 


-00  .-11  _  D0 

BoYo  c0  Y1  Do 


A?r2  +  bTM  ■  0  • 


(A2.4) 


(A2.7a) 


(A2.7b) 


(A2.7)  oan  be  solved  silmultaneously  and  the  results  are 


■ft 


■ft  -  W 


(A2.8) 


.0  0 

“  Vo 


B 


-1 


(A2.9) 


The  set  of  equations  obtained  from  equating  the  ooeffioients  of  6 
in  the  truncated  set,  (A2.2)  is 


bJyJ+c'S?  -0 

0r0  0  '1 


❖XrK1*  22+c;S2-° 


(A2.10) 


A0/1  1+B°y1+B_1r2  -0 

nrn-l  n'n  n  fn 


When  the  result  in  (A2.6)  is  substituted  into  (A2.10)  it  is  found  that 


Yn*°  for  n>3 


(A2.ll) 


and  the  equations  are  reduced  to  three  algebraic  equations: 


Vo^o  rl 


Al+Blrl+Bll7l+<:ilr2"  °1 


.0  1.--1  2 
*271+B2  7 2 


0. 


(A2.12) 


These  equations  oan  be  solved  simultaneously  and  the  results  are 


(A2.13) 


Y 


4 

2 


(Ajyj)/B 


—  1 
2 


Dico1“co1(Bivi+cilv2) 

aV1-  B°B_1 

Vo  B0B1 


-<bW,,c;1‘ 


(A2.14) 

(A2.13) 


Another  set  of  equations  oan  be  obtained  by  equating  the 
coefficients  of  62  and  the  set  is  shown  below. 


B°y2+C"1y3 
°0*(j  0  T1 


0 


A®y2+BBy2+B“lrj+C"lr2 


-0 


A2Yl4B2T2+Bal,y24C2lT3 


D 


2 

2 


A2  , 

nrn-l 


( A2.16) 


When  the  result  of  (A2.ll)  is  substituted  into  (A4.16),  a  generalized 
result  similar  to  (A2.ll)  is  obtained  for  r™  where  m-3. 


Tn"° 


for  n>4 


(A2.17) 


Then  the  set  (A2.16)  is  reduoed  to  just  four  equations  as  shown  below. 


B^y^+C”3* 
H0y(j  0  Y 


(A2.18) 


A2Y1+B2Y2+B2Y2+C1Y2^) 


,0  2  ,Q-1  3  A 
A3y2+B3  y3"° 


c 

2 


When  these  equations  are  solved  simultaneously,  the  results  are; 


1 


j)2_/ ,0  2  _0  2^,-1  3 
°2  (A2y1  B2y2^2  y3 


(A2.19) 

(A2-.20) 


(A2.21) 


(A2.22) 


Similarly  from  equating  the  coefficients  of  8  and  applying  the  result  in 
(A2.17) ,  a  general  result  for  yA. 

yJ-0  for  n>5  (A2.23) 

n 


The  set  of  simultaneous  equations  obtained  is 


❖taM  - 0 


*0^2 .«0  3  g—l  4^,-1  4  * 
AlWl  Yl4Cl  y2^ 


2  1  2r2  2  2  2  r3 


A3t2+B3y3+B31y3+C31y4"D 


3 

3 


a^tJ+b’VJ  -  0 


(A2.24) 


These  aquations  when  solved  simultaneously  give 


(A2.2S) 


D3"(A3T2+B3Y3+C31t4) 


B 


-1 


(A2.26) 


-fA°v3+B0v3+C"1v4) 
4  (A2yl  B2Y2^2  V 


B 


-1 


(A2.27) 


(A2.28) 

(A2.29) 


The  same  procedure  can  be  repeated  for  obtaining  higher  order  y®. 
In  general,  by  equating  the  coefficients  of  5,  it  will  be  found  that 

ym«0  for  m>n  (A2.30) 

n 

and  the  set  of  equations  that  is  left  to  be  solved  will  always  be  found 
to  contain  m+2  equations  with  an  equal  number  of  unknowns  as  shown 


(A2.31) . 


Vo 


+  c 


-1  ra+1 
0  Y1 


0 


2'1  2'2  2  '1  2  f3 


*SCi*Bft*B 


-1  n+l 
m  Ym 


+C 


m+1 

ra+1 


Vv>.- 


(A2.31) 


9.3  APPENDIX  3:-  ANALYTICAL  METHOD 

If  the  trunoated  set  of  equations  in  (2.10)  Is  assumed  to  have  four 
equations (i.e.  n«3),  they  are 


BoWi"Do  (A3-1) 

*1Y0+Blrl+C1Y2  "  °1  (A3. 2) 
A2y1+B2y2+C2y3  "  °2  (A3. 3) 
*3y2+B3y3  "  B3  (A3. 4) 


where  the  coefficients  A,  B,  C  and  D  are  defined  in  section  2.  These  set 
of  equations  can  be  solved  simultaneously  and  the  results  are; 


(M.5) 

Yj  »  ^B2B3“C2A3^D1_A1y0,“C1^B3D2“C2D3^ 

(B2B3“C2A3)BrB3A2  (A3.fi) 

y2  "  B2D3"C2D3_B3A2y1 

B2B3-C2A3  (A3. 7) 

t3  -VVi 

B3  (A3. A) 

The  coefficients  A  ,  B  .  C  and  D  by  A0,  B°+B“1/8,  C-1/5  and  DnSn 

nnn  n  nnn  n  n 

respectively.  Those  terms  with  power  of  S  greater  than  -2  are  assumed 


6 


T 0  .{‘■3<DSBrlBilB3l’+‘'JtD0(BlBila3lt8IlB2lBrAJO^Bi1,-DlC0lB2lB3,,j 

K3‘B^VBI1-*JCjVB31,+*‘JtBS(B!BIlB31*BIlB2B31*BIlB2lB3,.1 


*1  . Vl*?YOt02lA3-BX1-BrBrBIlB31‘'11  *B2lBIlD!r1} 


* 


•■3,BlVh‘,t,'i(BXlB3l*BIlBSBJ,*BrB^BrD7lcI1‘3> 


.8'l(BjB'IBj*Bj1BjBj*BjBjB"1-Bjc"l»j)p 


2  - 


(DrA3Y2)/(B3+B31) 


(A3. 10) 


(A3. 11) 
(A3. 12) 


If  a  new  variable,  y.  Is  defined  such  that 


y»  (?/o)1/2  (A3. 13) 

then  the  coefficients  can  be  rewritten  as 

aJ  -  C"1  «  1  -(2y2+l)  sinh~1y  (A3. 14) 

Ay2(y2+2) 

A?  -  C71  »  1  -(2y2+3)  sinh“1(l+y2)1/2  (A3. 15) 

IT,  2+l)(y2+2) 

k°3  -  C"1  -  1  -(2y2+3)  sinh~1(y2->-2)1/2  (A3. 16) 

/( yZ+2)(yZ+3) 

BB  -  B^1  -  2  sinh”1y  (A3. 17) 

B?  -  Bl1  «  2  slnh_1(y2+l)1/2 


(A3. 18) 
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bJ  -  B"1  -  2  sinh-1(y2+2)1/2 


B^  -  2  slnh_1(y2+3)1/2 


D"-  -<r(y2+n)e“  1 


Taylor  expansions  oan  be  used  to  expand  the  Inverse  hyperbolic 
functions.  In  this  case  the  expansions  are  made  about  y-O. 


sinh“1(y2+n)1/2  -  sinh"1(n)1/2  + 


(A3. 19) 


(A3. 20) 


(A3. 21) 


■v  -  jC2JL2n±U,  +... 
Vr?+l)  $24  (n'+l)^2) 


(A3. 22) 


By  using  the  expansions,  the  A0,  B®,  B-’1  and  C-1  oan  be  written  as  a 

n  n  n  n 

series  In  power  of  y  as  shown. 


Ai  *  Gnl  “  4y2(~1/3  +  y2/5  ♦•••• 


*2  "  Cll  *  -0.5942y2  +... 

A3  "  C22  *  -1*3398  -0.4174y2  +... 


bJJ  -  B”1-  2ytl  -y2/«  +3y4/40  +... 
bJ  -  Bj1-  2tsinh~1(l)  +y2/(2 SI)  +.,, 


bJ  -  B"1-  2tslnh_1/T  +y2/(2/T)  +... 


bJ  »  2tslnh”1/T  +y2/10  +...  ] 


(A3. 23) 


(A3. 24) 


(A3. 25) 


(A3. 26) 


(A3. 27) 


(A3. 28) 


(A3. 29) 


These  expansions  can  be  substituted  Into  the  expressions  for  gQ,  and 
and  the  final  approximations  are  shown  In  equations  (4.7),  (4.8)  and 

(4.9)  of  section  4. 


When  deriving  the  analytical  formula  the  formulas  derived  from 


•-V-V-V-I.'l 

'  *.  - 

> « v  *.* 


-V.V-; 


viy.vV-d 


■V vv%>>, 

,  *.  *.  %  v1 


i*.  •*  •' %  *.***•.  «_“***•**•  «  *■*«  *  «  •  »  m  *“K«J 


y  Vf-sA 


Perturbation  Method,  the  coefficients  are  expanded  in  terms  of  y  using 
equations  (A3.22-A3.27) .  In  this  ease  the  different  terms  in  rn 
expansions  used  are  given  in  equations  (3,13-3.20). 


-102 


9.4  APPENDIX  4:-  EXPANSIONS  FOR  HIGH  ELECTRIC  FIELDS 

Fop  a  parabolic  band  fl«£and  if  there  is  no  p-wave  scattering  ,  e*0 
and  a-1  in  equations  (A1.10)  and  CA1.11>.  These  equations  then  reduce 
to  equation  (5.3.7)  and  (5.3.8)  in  Conwell's  book  T251. 


(do  -  «E„H  { 

'dt  '  c  syo  7  (in*h(dx) 

+H  (x— 1 )  t  f  (x-1 )  -e^f  n  (x )  ]$J(x- 1  )1 

0  0  J  (A4.1) 

dk^i\  -  eE„cos0N/e',g  (x+1)  /fx+1)  -5]  -g(x)^x) 

**  -f -  \  V<~)  L  277x7x717)  J  r 


m 


asy 


+H  (x-l)  {a  (x-1) _/(x-l)[(2x~1)^(x“1)-j-g(x)eq(x-l)|^ 


f(  x  )L2/x(x-t) 


(A4.2) 


where  the  terns  have  been  replaced  by  ^(x).  The  rate  of 

change  of  the  distribution  function  due  to  the  applied  electric  field, 
(A1.12) ,  is  given  by 


/df(J|  -  -eSf  2  d(x3/2g)  +  hkcosftif  1  (A4.3) 

\dt  J?  I3h/x  dx  ra*  dx  ' 


When  (A4.1)  is  equated  with  the  symmetric  part  of  (A4.3),  and  (A4.2) 
equated  with  the  asymmetric  part  of  (A4.3),  the  resulting  equations  are 
shown  below. 

d£xfj(x)  1  „  n^i[f0(x)-e,Tf0(x+l)]  ^Xx)-H(x-l)[f0(x-l)-e<,f0(x)]^(x-l)} 

( A4 .4) 


xdf^(x) 

dx 


1  (fj 

3X\  1 


(x)^(x)-e'7f1(x+l) 


(2x+l)0(x) 

27xTx+TT 


(A4.51 


where  fj(x)  -kg  and  E«2E/(3Eq).  It  la  assumed  that  the  electron 
distribution  funotlons,  tf 
power  series  In  1/E  where 


distribution  funotlons,  f^Cx)  and  fj(x)  can  be  expanded  into 


f0(x)  "  E<N/K)J  f£(x) 
fl(x)  -  Z(N/R)J  fi{x) 


(k4.fi) 

(A4.7) 


Therefore  (A4.4)  and  (A4.5)  become 


d_[x£<l/K>JfJ<x>]-s/K  [£a/K)Vj(x)^"£(l/K)Vj(xM>]<£(x> 

-H(x-1)[  £(l/E)Jf^(xHOVT£(l/KV1fJj(x)|0(x-l) 
J  j  (A4.R) 


ax  1  1  1  ■ 


+H(x-1)[  e^(x-l)^(l/K)  jfJ(x) 

(JM.„ 


where  the  superscripts  of  (1/K)  is  its  power  while  the  superscripts  on 
f0(x)  and  f^(x)  are  Just  there  to  Indicate  the  order  of  each  terms  in 


the  expansions.  The  coefficients  of  the  different  are  equated  and  the 
equations  that  are  found  for  the  first  few  power  of  (1/K)  are  shown 
below. 


For  (1/K) 


x  <LLfo(x)J  "  0 

dx  ° 


d[xf,(x)l  -  0 
dx 


(A4.10) 

(A4.ll) 


For  (1/K) 


dCxf,  (x)I 
dx  1 


hj[fJ(x)-e,Tf2<x+l) 


O'" 

-H  <x-l  )0<x-l ) ff2 (x-1 )  -e^f  2  (x  )1  j 


(A4.12) 


’(x+1) 


xd(fA(x)1  -  N;$(x)f!;  (x)-e*  K2x+l)&x)  -i]  f( 
dx  0  3 1  |_7x7x?n  J 

+H (x-1 )  je^x-l ) f J (x ) -  |j2x-l)4(x-l) -tj  f J (x-ljj 


( A4.13) 


For  (1/K)‘ 


dixfj(x)]  -  M|(fJ(x)-e'TfJ(x+l))(^x)-H(x-l)^Cx-l)(fJ(x-l)-«fffJ(x|j} 


14) 


15) 


For  (1/K) 


d(xfj<x)1-lf{|^(x)-«"f*(x*l)]0(x)  -  H(x-l)#x-l)[fJ(x-l)-e°fJ(xft 

T A4 • 

xdifj(x),.  sj  fr,)rj<x)-.’[(g%jffij|  fj(x-l) 

+H(x-l){«'T0(x-l)f?(x)-  j^>(n?)  -ll  fj(x-l)l! 

1  r  1  J  1  f  |  (A4. 


16) 


17) 


From  (A4.10),  for  the  aquation  to  be  true  for  all  x,  fQ  has  to  be  a 
constant  and  Is  assumed  to  be 

fj|<x)  -  FO)  (A4.1*) 

From  (A4.ll),  it  is  found  after  integration  that  xf^(x)  is  a 
constant.  By  taking  the  limit  as  x  approaches  zero,  it  can  be  seen  that 
the  Integration  constant  has  to  be  zero  or  else  f^(x)  is  infinite  at 
x»0. 

m  fj(x)  -  0  (A4.19) 


This  result  is  substituted  into  (A4.14)  and  it  is  reduced  to 


xd  f*(x)  -  0 
dx"q 


(A4.20) 


For  this  equation  to  be  valid  for  all  x,  f*  has  to  be  a  constant.  For 
fQ(x)  to  match  the  result  of  Stratton  [J],  and  Conwell  and  Vaasell  [dl 
at  very  large  B  and  x,  f®(*)«^,(0).  This  last  result  implies  that  all 
for  all  m>0.  Therefore  the  oonstant  fj(*)  must  be  zero. 


fJ(x)-0 


(A4.21) 


Equation  (A4.18)  can  be  substituted  into  (A4.12)  and  integrated 
for  an  expression  of  f*(x>  in  terms  of  F(0)  and  x 


f*(x)-F(0)  (l-eq)|2 J’ainhT^/x  dx  -2 J" sinh-1/x-l  dx| 


(A4.22) 


f*(x)  is  zero  at  the  lower  limit,  x«0.  In  the  second  integral  the  lower 
limit  is  the  result  of  the  step  funotion.  The  method  of  substitution  is 
used  to  Integrate  the  inverse  hyperbolic  sine  function. 


1  rfv 

Let  y«/x,  therefore  dy-j 


Then  / sinh“Vx  dx 


I 

2y  / sinh_1y  dy 


■  °<x+l/2)sinh"Vx  -  —SEP 


(A4.23) 


Equation  (A4.22)  after  the  integration  is  given  by, 
f*(x)  -  F(0)  (l^q)|(2x+l)alnh~Vxr-/x(x+i) 

-H(x-l)[(2x-l)sinh”1yx-Wxtx-l)j  } 


(A4.24) 


This  equation  can  be  substituted  into  (A4.15)  and  the  resulting  equation 


oan  be  integrated 


+fJ(0)  (A4.25) 


2  2 
where  Fg  is  a  constant  of  integration  that  will  reduce  fg(x)  to  zero  at 

x—».  If  f^(«)“0,  then 


F*  (x) — 1 / dx/2f« (x) sinh 


Til)1 


l/x  -  e^f.  (x-»l>  fex+fl  slnh~Vx  -ll 

1  ["Tx7x+1)  J 


-H(x-l) !  2ea  fhxJsinh”1/^  - 


f|(x-l> 


:-0  ainti/^T^I  -l|  i 

A(x-i)  J) 

(A4.24V 


From  equation  (A4.14),  f^(x)  can  be  written  in  the  integral  form 


f*(x>« 


X 

1,  /dx{[f*(x>  “ 

•  -H(x-l)  ^fg(x-l)  -  e^f^x)]^-!)! 


(A4.2A) 


This  equation  can  be  simplified  by  Integration  by  parts.  Consider  the 
first  part  with  #x)  replaced  by  2sinh-1/x, 


I  -  Jdx  |fj(x>  -  •  fjj(x+l)j  2slnh 


Let  dr  -  2sinh-1/x  dx 
and  a  ■  tfJ(x)-e*fj|<x+l)] 
therefore. 


I  •  £f*(x)-e*f*(x+l)J£  (2x+l)sinh”Vx  -  /x(x+i)J 

-  /u2x+l)alnh“l/x  * 

•l  (A4.29) 

where  the  result  from  (A4.23)  has  been  used  for  the  integration  of  the 

inverse  hyperbolic  sine  functions.  Therefore  (A4.28)  can  be  written  as 


fj(x)-l/x  jK(x)“*‘,fo<x+1>]  [<2x+l>sinh“1/x-/x(x+l)] 

-H(x-l)  [f5(x-l)^fJ(x)][(2x+l)ainh“1/3F^-/xTxTT^ 
- Jix  [(2x+l)sinh~Vx-/xtx+i)J  Jctrjcx>— dfj(x+l)j 
♦ Jix  J(2x-l)sinh~1/x-l-/x(x-l)|  |dT^ (x— Cx )|  | 


(A4.30) 


>2a<Tg(g,)pc,2-Mc2)(k,2-Mc2  -l)  ln|k'+k}  -k'/k^+k2 

k  {.  2k  A  2k 'k  '  Ik’-kl  \  kk’ 

J2g(g)  |/k'2+k2  -2Un|k^+k|  -2l  dfl 

U  kk’  /  Ik'-kl  J  d£ 


[k'2+k2ln 

k'+k 

-k'" 

-e^gdJln 

k'-fkj 

[  2k 

k'-k 

k'-kl 

*g(D  |k*2+k2fk'2+k2  -iiln|k»+k{  -k,/k,2-»-k2\-2"| 

!c  L  kk'  \  2kkr  /  Ik'-kl  \  kk'  /  j 

*«"g(t)[7k'2+k2  -2\lnlk'-k|  -2]  do  ) 

H  kk'  /  ik'-kl  I  dt  [ 


In.  the  two  equations  above  terms  that  are  assoolated  with  p-wave 

2 

scattering  are  terms  with  c  .  If  there  is  no  p-wave  scattering,  then 
o*0,  and  the  equations  will  reduce  to  (A4.1)  and  (A4.2). 

It  is  assumed  that  the  non-parabollolty  of  the  band  is  of  the  form 

fl-fc+nfc2  (AS. 3) 

therefore  dD  _ 

d V  UUt 

The  next  step  is  to  express  the  k  and  k  *  in  terms  of  t and  tiu  by  using 
the  relation 

T>*k2  -£M£2  (AS. 4) 

2  a* 


k'*+k* 


k'k 


7<i+ee*WMle+a(i+TleJ*J 


l+2*La 

(£+o<?)  J i  ♦2dhwt+g(ii%i^-»'n« 
t+«tJ 


-  2+ 


1  r 2ahw£-*q(hw)1-»ft<»,|»-f  .  . 

7  L  j 


For  £>>hw  this  oan  be  approximated  by 


k*»+k*|  -2  +  i  r2ohMg»hw|2 

k'k  j  4  I  e+«t* 

f'-C+hw  J 


(AJ.J) 


Similarly  for  the  oase  wnen  £'-£-he,  the  same  procedure  gives 


k*»+k»{ 
kk'  I 


-2  ♦  1  /-2ahw£--h»\* 


(AS. 4) 


_  u  •  Air 

The  term  «an  *>•  written  in  the  form  of  inverse  hyperbolic  sine 


functions  and  are  given  by 


In  k'+k 


In  k*+k  -  2sl 


i  k»+k 
k'-k 


£>£*&•» 

for  the  case  when  £•  -£+h« 


yk'*-k» 

t >#*» 


KKi.TinfdB.Tii 


for  £>>h*« 


Vha  1+533 


(A5.9) 


Similarly  for  the  ease  when  £'*£rh«, 


k* 

7ki-k'* 

f-l-h- 


_  rg  (1+el)  .  r 

•.2l«<l+2oc)  ”  I 


(A5.10) 


when  equation  (A5.4)-(A5.10)  are  substituted  into  (A5.1),  and  only  the 
terms  due  to  p-wave  scattering  are  considered,  the  equation  can  be 
written  as 

df0(l)  «eWEqc»j  |(l+2a£+2oh«)  (f0(*‘*l»)e<r-f0(C))| 

r  ->•} 

+tl+2a£-2afte)  (fQ  «-<»*> -e"f0<«>)j 

/lPahn  +h'0l».J_^-ir  £  (l+a£)  .11,Z  J 

\vrz+*v  ]  sinh  IhSd^)-1!  ’2f 

di5.ii) 

ut  th">  tlM  two  inverse  hyperbolic  sine  terms  can  be 

written  as  slnh~Vx  and  slnh~*/x-l .  Since  the  case  considered  here  is 
2x*  ix* 

for t  >>tk*,  x  is  quite  large  compared  to  unity.  However  for  large  x, 
sinh“Vx  is  rather  small  oompared  to  2. 


w  >  +  ■  «  »  * 

h»(l+2ac>  1 


(A5.ll) 


df  (C)  — 2eE.no*  j  (l+2a£+2ah*»)  (f. (£+h«)eff-f.(e)) 

dtu  hk  '  u  u 

+(l+2at-2ah»)  (f0(£-'h»)-e<,f0(l))|  (AS. 12) 

Taylor  expansions  are  aade  about  for  f^d-Fh*)  and  if  only  the  first 
three  teras  are  retained,  the  equation  oan  be  simplified  to 

2eE0e»  |  4Oh«f0(t)  +  hi^df 0  ( 1)  [l+2o£+  ( 2N+1 )  2aftiJ| 
hk  ( 

+(h»»)*^f0(f )  [(l+2a£)(2N+l)+2aT)»]|  (A5.13) 

When  p-wave  scattering  is  Included,  the  Boltzmann  equation 
describing  describing  the  symmetrical  part  of  the  distribution  function 
in  the  lower  valley,  assuming  no  intervalley  scattering  can  be  written 
as  f  61 

***■  f-T'-,11/2j  rrt*W/2  x  f„  i  ♦  /df„\  *  .  „ 

3m*(l+2al)U+o£2l  d*U.  hw  I  W«t  I  \ dt  '  c-sw'dt  '  c-pw 

(AS. 14) 

where  the  non-parabolic ity  in  the  energy  band  of  the  form  in  (AS. 3)  has 
been  used  for  the  field  term.  The  collision  term  has  also  been  written 
as  the  sum  of  a  s-wave  scattering  and  a  p-wave  scattering  terms  as  given 
in  (7.3).  x  is  the  polar  optioal  relaxation  time  for  £.>rhw  as  discussed 
in  section  6.  For  the  non-parabolieity  considered  in  this  section,  it 
is  given  by  (51 

T  _  r2m*h«11/2  e*  -  l  (l+*V)in  (AS. IS) 

l  eEn  I  9*  +  1  J (hw)  (l+2ot) 


The  collision  tern  dfft(e)  is  identical  to  the  result  when  no  p-wave 

dt 

scattering  is  considered.  For  the  type  of  non-parabolic ity  considered 


here  it  oan  be  written  as  [29] 


dfA 

3t°o 


eE(htt) 


—  + 
hk{l+2afc)di(lh« 


<W+l/2>df 


(l+2ot)ln 


fi 

Lhtt  l+2«£jj 


(AS.lfi) 


Equation  (A5.14)  can,  therefore,  be  written  as 


2c»E» 

3eE0n* (l+2at) 


r  2b«  r  ^/e*  -  i\  dr  ( t+o£_>  d  f  0 1 

U+oi2  I  \e'T  +  l/dtl.(l+2aOdt  I 


+  eE.fhtt) 


Dk(l+2«£) 


*  d Ljrffl  +(N+l/2)dffl-J(l+2«t)*lnr4_(t+n£ 

.)  d  gjLh'i  dll  Lhw(l+2a£). 


„  *2eE-0*r4otJ«fft  +  •h«dfA(l+2ol+(2N+l)*ah») 


d’f  ( 


(l+2af  ) (2N+1) +2ah«) 


T  i? 


(A5.17) 


It  is  assumed  that  the  solution  of  A5.17  can  be  written  as  the  sura  of  a 
dominant  term  f®(x)  due  to  s-wave  scattering  and  a  small  perturbed  term 
«fj(x)  caused  by  p-wave  scattering. 


#0  x  A 

fo  *  ,ro 


(AS. 18) 


where  s  is  a  smallness  parameter  to  be  set  to  one  later.  When  this  sum 
is  substituted  into  AS. 17,  the  terms  in  the  equation  can  be  reordered  in 
terms  of  and  «*.  The  coefficients  of  each  power  of  «  can  be  set  to 
zero  and  two  equations  are  obtained. 


.»  *  •  *  •  *  •  *  •  •  •  *j  *,»  •  •  *  m  *»  •  -  *  »  ’/  '/  *  *  *  •  *_*  •  %  V  *j  ^  ^  '•*  V  V  *  ■  *  *  V*  ' 


2e»E* 


_ pm*  r  !•  -A<*  rtg-mc1)*  df,,] 

3«E0a»(l+2ot)  U+ol2  I  xe^+l  /dU(l+2«t)2  d£  J 

+  !V^I  dK  +(M+l/2)  1£1(  i+2o£)  *  in  CL  1  j 

-hk(l+2«£)d^|th**  d£.'  fow(H>2«t)  J) 

(A5.19) 

2e»B*  f  2m*  ]1/2  eg  -1  d  [,(t+at>>> 

3«B0m*(l+2«£)U+a€.2 1  e*  +1  dJ.(l*2«!)i  it  I 


Wc(l+2ol)d€(U«  dl  l  bl«(l+2o£jJ) 


2eE„o* 


h«!£o  (l+2a£.) 

d  e 


(he)  *  — *0 ( l+2o£)  (2N+1)  1 

2  d£2  J 


(A5.20) 


The  first  of  these  equations  is  the  one  that  is  obtained  by  equating  the 
eoeffioients  of  the  first  Legendre  polynomial  PQ  to  zero  in  the 
Boltzmann  equation  with  p-wave  scattering  neglected.  This  equation  has 
been  solved  by  Conwell  et  al  [51  for  an  arbitrary  energy  band.  The 
solution  that  they  derived  is  shown  below. 


F(0)  exp- 


fe1 


4  lnMn/y)  dx  1 
+  d(^04  ln(4q/q 


(AS. 21) 


where  x»5f he  as  defined  in  seotlon  5,  b  and  d  are  defined  by  equation 
(5.14)  and  (5.15)  in  section  5.  n  is  the  energy  function  D  normalised 
by  the  he.  For  the  type  of  non-parabollolty  considered  here, 
q»(€+a£*  )/(!)*)  end  the  prime  in  n  Indicates  differentiation  with  respect 
to  x.  F(0)  is  a  oonstant  of  Integration. 
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To  see  how  much  p-wave  scattering  affects  the  electron  distribution 
function,  equation  (A5.20)  has  to  be  solved  for  f*(x).  There  are  three 
terns  within  the  bracket  that  are  multiplied  by  the  factor  c*.  The  last 
of  these  three  terms  is  a  product  of  the  second  derivative  of  f®(x)  and 
fhw)s,  whloh  is  small  oompared  to  the  other  two  terms  and  is  dropped. 

2e*B»  f  2m*  11/2  eg-l  d  r(t+n6.»)»  df jh 
3eE0m*(l+2«Ol£+otZJ  e^+l  dil(l+2o€)2  d&  I 

hk(l+2a£)d  E-jlnu  it  '  Lh«(l+2otjJ) 


ft«lTo(l+2o£) 

it 


0  . 


(AS. 22) 


This  equation  can  be  written  in  a  more  compact  form. 


«» 

3E_ 


( fi+ 


e*-l  d  r  g+«£»  dfll  %,d 
e^l  dtl.(l+2ot)»dt°i  (n“  d£|hS 


(N+l/2)df; 


:) 


~4  (  E+*  g*  )  1 


-0  (i+2«0*ln  S^^j 


-2ft«c»  df(l+2o€)*f[J]  -  0 


d? 


(AS. 23) 


In  terms  of  the  variables  x  and  n  (AS. 23)  can  be  written  as 


-0 

dx 


(A5.24) 


This  equation  can  be  integrated  to  give 


dfj+  f2n>4ln(4n/n’) 
dx  b^  +drj'  ln(4W) 


,  91 


bn  +dn'  ln(4n/n’> 


£ 


*d[n'*fg  1  dx 


(A5.25) 


In  this  equation,  o*  is  also  a  function  of  t,  and  .has  to  be  written  in 
terms  of  x  so  that  the  right  hand  side  can  be  Integrated,  o*  can  be 
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written  as  Ca3! 


c  4 


_  i 


JL 


1/2 


whs re  gQ  is  given  by 


to  -MS? 

2m* 


(A5.26) 


(A5.27) 


••■•v-v 


where  m^  is  the  mass  of  electron  in  vaouum.  Since  is  related  to  a* 
by  m_«m*  (1+ (£„/&.)) ,  e*  can  be  written  as 

o  p  g 


o*  4 


n»k*  4g 
1+  ■  p 


(A5.29) 


Since  h*k*/(2ra*)  is  Just  £+«€.*  for  the  type  of  non-parabolicity 
considered  here,  c*  can  be  rewritten  in  the  form 


o>»‘  ;i- 


1+  (W1 

WV«p 


1/2 


(A5.29) 


For  most  of  the  compounds  considered  here  t  is  much  larger  than  Z  . 

P  g 

Therefore  the  <£  +£  )/£  can  be  approximated  by  unity.  The  factor 

r  o  r 

£*/t£^(£g+£p)]  is  Just  a  t<?T.  Therefore  c*  oan  be  approximated  by 


c*«J(l-l/i,'> 


(AS. 30) 


This  expression  is  substituted  into  the  right  hand  side  of  (A5.25)  and 
integration  by  parts  is  used  to  integrate  the  left  side. 
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dx  bn  +dn#  ln(4n/n')  bn  +dn'  ln(4n/n 


Pr on  (AS. 21)  it  is  found  that 


(A5.31) 


d_Jln(F(0)/f”)> 

dx 


f  n>ln( 
IbV+dV 


’lnMn/n*) 
+<1ti  'ln(4ri/n  1 


(AS. 32) 


Therefore, 


diD.folL1"(P<0,/fO)1. 

dx  dx 


bn  +dn*  ln(4n/n 


f?-rS",-2'‘,wr  **) 
•) '  ' 


(AS. 33) 


It  is  possible  to  transform  the  right  hand  side  of  the  above  equation 
into  the  form  of  an  exact  integral  by  the  integrating  factor,  F(0)/f®. 
With  the  Integrating  factor,  (AS. 33)  can  be  integrated  from  a  lower 
limit,  x,  to  x.  This  lower  limit  is  chosen  such  that  fj/f®  vanishes  at 
x*A.  The  lower  limit  is  found  to  be  one.  After  the  integration,  F(0) 
cancels  out  from  both  side  of  the  equation,  and  the  result  of  fQ(x) 
expressed  as  a  ratio  of  fj  is 


-» -  ff-r1-* - kv-’tJfjV 

f0  Jlbq  +dq’4ln(4i,/T,»)  '  fji 


0  *« 


(AS.  34) 


A 

For  x  smaller  than  one,  bn*  is  much  smaller  than  dV  ln(4n/n’). 
Therefore  for  the  region  0<x<l,  f®  can  be  approximated  by 


fj-  P(0)  exp 


-/l/d  dx 


-  m  d( 1-exp H/d)). 


This  allows  (AS. 34)  to  be  written  as 


W„N\% 
*  s'  V*  % 
,*  Vv 

.•v  *  *N 


/  /V.V 


*  *-n '-^(p  (o«  i-e"1'*1 ) Jdx)|Jdx 

(A5.3S) 

For  large  eleotrio  field  and  x,  tha  asysnetrloal  part  of  the  electron 
distribution  function  In  teras  of  fQ  Is  [5,6] 

ft  ,  -hegT  dffl  (AS. 36) 

1  a»hwi»  dx 

where  x  Is  given  In  (A5.16).  In  terns  of  f^  and  f*,  f^  is  given  by 


kg  -  j*>cgT  r 

«*  ♦  . 

(«S » *  i£  f  Pdx\i 

m*Ti**r| 

dx 

'  dx  ■'*  '  1 

(A5.37) 
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